On Features and Nongaussianity from Inflationary Particle Production 
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Interactions between the inflaton and any additional fields can lead to isolated bursts of particle 
production during inflation (for example from parametric resonance or a phase transition). Inflation- 
ary particle production leaves localized features in the spectrum and bispectrum of the observable 
cosmological fluctuations, via the Infra-Red (IR) cascading mechanism. We focus on a simple proto- 
type interaction g 2 (cj> — 4>o) 2 X 2 between the inflaton, <j>, and iso- inflaton, \\ extending previous work 
on this model in two directions. First, we quantify the magnitude of the produced nongaussianity by 
extracting the moments of the probability distribution function from lattice field theory simulations. 
We argue that the bispectrum feature from particle production might be observable for reasonable 
values of the coupling, g 2 . Second, we develop a detailed analytical theory of particle production 
and IR cascading during inflation, which is in excellent agreement with numerical simulations. Our 
formalism improves significantly on previous approaches by consistently incorporating both the ex- 
pansion of the universe and also metric perturbations. We use this new formalism to estimate the 
shape of the bispectrum from particle production, showing this to be distinguishable from other 
mechanisms that predict large nongaussianity. 



PACS numbers: 11.25.Wx, 98.80.Cq 



I. INTRODUCTION 

The inflationary paradigm has become a cornerstone 
of modern cosmology. As measurements of the Cosmic 
Microwave Background (CMB) radiation grow increas- 
ingly precise, it has become topical to look beyond the 
simplest single-field, slow-roll inflationary scenario. In 
particular, it is interesting to determine the extent to 
which non-minimal signatures, such as features in the 
primordial power spectrum or observable nongaussian- 
ities, can be accommodated by microscopically sensible 
inflation models. Efforts in this direction are valuable be- 
cause they allow us to test our theoretical prejudices and 
provide observers with well-motivated templates for de- 
partures from the standard scenario. Finally, a detection 
of some non-minimal features might open a rare obser- 
vational window into fundamental particle physics at ex- 
tremely high energy scales. In this work, we will consider 
a very simple and well-motivated class of models which 
predict novel observable signatures in the spectrum and 
bispectrum of the primordial curvature fluctuations. 

In a variety of inflation models, the motion of the in- 
flaton can trigger the production of some non-inflaton 
(iso-curvature) particle during inflation. Models of this 
type have attracted considerable interest recently; exam- 
ples have been studied where particle production occurs 
via parametric resonance P-Q, as a result of a phase 
transition 

EMI, 

or otherwise [l8|. Such constructions 
are novel for a variety of reasons: 

1. The produced iso- inflaton particles may rescatter 
off the slow-roll condensate and generate a signifi- 
cant contribution to the primordial curvature fluc- 
tuations through the process of Infra- Red (IR) cas- 
cading pj. This provides a new mechanism for gen- 



crating cosmological perturbations that is qualita- 
tively different from the standard mechanism p^ |. 
the curvaton [13, Hl| or modulated fluctuations 

Hi. 

2. Particle production and IR cascading leads to a 
variety of novel observational signatures, including 
features in the primordial power spectrum and also 
nongaussianities [l], 0. 

3. Particle production arises naturally in a number 
of microscopically realistic models of inflation, in- 
cluding examples from string theory [3| and su- 
persymmetric (SUSY) field theory [27| • In partic- 
ular, inflationary particle production is a generic 
feature of open string inflation models 0], such as 
brane/axion monodromy (28l - [3l| . (See also [Hj].) 

4. Observable features in the primordial power spec- 
trum, generated by particle production and IR cas- 
cading, offer a novel example of the non-decoupling 
of high scale physics in the Cosmic Microwave 
Background (CMB) [111]. In the most interest- 
ing examples, the produced particles are extremely 
massive for (almost) the entire history of the uni- 
verse, however, their effect cannot be integrated out 
due to the non-adiabatic time dependence of the 
iso- inflaton mode functions during particle produc- 
tion. In Q particle production during large field in- 
flation was proposed as a possible probe of Planck- 
scale physics. 

5. The energetic cost of producing particles during 
inflation has a dissipative effect on the dynamics 
of the inflaton. Particle production may therefore 
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slow the motion of the infiaton, even on a steep po- 
tential. This gives rise to a new inflationary mech- 
anism, called trapped inflation 0, [33[ , which may 
circumvent some of the fine tuning problems asso- 
ciated with standard slow-roll inflation. See [3] for 
an explicit string theory realization of trapped in- 
flation and [lij for a generalization to higher dimen- 
sional moduli spaces and enhanced symmetry loci. 
The idea of using dissipative dynamics to slow the 
motion of the infiaton was exploited also for a very 
interesting mechanism (which pre-dates trapped in- 



flation) called warm inflation 
variant of natural inflation [35 
recently by Anber & Sorbo [18 1 



|27J. See also the 
that was proposed 



In this article we study the impact of isolated bursts 
of inflationary particle production on the observable pri- 
mordial curvature perturbations. In order to illustrate 
the basic physics, we focus on a very simple and general 
prototype model where the infiaton, 0, and iso-inflaton, 
X, fields interact via the coupling 
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On physical grounds, we expect that our results will 
generalize in a straightforward way to more complicated 
models, such as fermion iso-inflaton fields, gauged inter- 
actions and (perhaps) inflationary phase transitions. 

Scalar field interactions of the type ([T]) have also 
been studied recently in connection with non-equilibrium 
Quantum Field Theory (QFT) [37]], in particular with ap- 
plications to the theory of preheating after inflation |4l| - 
l46j and also moduli trapping [33|, [34[ at enhanced sym- 
metry points. Although our focus is on particle produc- 
tion during inflation (as opposed to during preheating, 
after inflation) some of our results nevertheless have im- 
plications for preheating, moduli trapping and also non- 
equilibrium QFT more generally. For example, in [I| 
analytical and numerical studies of rescattering and IR 
cascading during inflation made it possible to observe, for 
the first time, the dynamical approach to the turbulent 
scaling regime that was discovered in [47|, |48( . 

Let us now discuss briefly the physics of the model ([!]) . 
At the moment when <f> = <po (which we assume occurs 
during the observable range of e-foldings of inflation) the 
X particles become instantaneously massless and are pro- 
duced by quantum effects. This burst of particle produc- 
tion drains energy from the condensate <fr(t), temporarily 
slowing the motion of the infiaton background and vio- 
lating slow roll. Shortly after this moment the \ particles 
become extremely non-relativistic, so that their number 
density dilutes as a -3 , and eventually the infiaton settles 
back onto the slow roll trajectory. 

The dominant effect of particle production on the ob- 
servable spectrum of curvature fluctuations arises be- 
cause the produced, massive \ particles can rescatter 
off the condensate to generate bremsstrahlung radiation 
of long-wavelength 8<f> fluctuations via diagrams such as 



Fig.[TJ Multiple such rescatterings lead to a rapid cascade 
of power into the IR. The infiaton modes generated by 
this IR cascading freeze once their wavelength crosses the 
horizon and lead to a bump- like feature in the primordial 
power spectrum. This bump-like feature is accompanied 
by a localized, uncorrelated nongaussian feature in the 
bispectrum 
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FIG. 1: rescattering diagram. 

In this paper we extend previous work [f], 0] on the 
model (JTJ) . First, we re- visit the problem of quantifying 
the magnitude of the produced nongaussianity. Using 
lattice field theory simulations we compute numerically 
the skewness and kurtosis of the Probability Distribu- 
tion Function (PDF) of the primordial curvature fluctu- 
ations. By comparison to the more familiar local model 
of nongaussianity, we argue that the bispectrum associ- 
ated with this mechanism may be observable in future 
missions. 

Next, we provide a detailed analytical theory of the 
quantum production of \ particles and the subsequent 
rescattering off the slow- roll condensate for the model ([IJ . 
This new formalism improves significantly upon previous 
efforts [1] by consistently incorporating both the expan- 
sion of the universe and also metric perturbations. We 
test our approach by comparison to fully nonlinear lat- 
tice field theory simulations, finding excellent agreement. 
We also use our formalism to estimate the shape of the 
bispectrum. 

The outline of this paper is as follows. In section \H\ we 
review the key results of P, 0|, describing heuristically 
the underlying mechanism of IR cascading and the resul- 
tant observational signatures. In section lrnl we character- 
ize the size of the nongaussianity associated with particle 
production and IR cascading, relying primarily on lattice 
field theory simulations. In section HVl we provide an an- 
alytical theory of inflationary particle production and IR 
cascading in the model ([1}, neglecting metric perturba- 
tions. Using this new formalism we estimate the shape 
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of the bispectrum in the model (JXJ) . In section [V] we 
reconsider our analytical approach, showing how metric 
perturbations can be consistently incorporated and, fur- 
ther, we demonstrate explicitly that their inclusion does 
not significantly alter the results of section IIVI Finally, 
in section IVH we conclude. 



where v = 0(0) and we have arbitrarily set the origin 
of time so that t = corresponds to the moment when 
(j> = 0o ■ The interaction ((T|) induces an effective (time 
varying) mass for the \ particles of the form 



(6) 



II. OVERVIEW OF THE MECHANISM 

In this section we provide a brief overview of the dy- 
namics of particle production and IR cascading in the 
model ([T]) and we also summarize the key observational 
signatures. This section is largely review of [H, 0], the 
reader already familiar with those works may wish to 
skip ahead to the next section. 

We consider the following model 
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where R is the Ricci curvature constructed from the met- 
ric g^, <fi is the inflaton field and \ is the iso-inflaton. 
As usual, we assume a flat FRW space-time with scale 
factor a(t) 



ds 2 



g^dx^dx" = -dt 2 + a 2 (t)dx 2 



(3) 



and employ the reduced Planck mass M p = 
(8ttG n )^^ 2 ^ 2.43 x 10 18 GeV. We leave the poten- 
tial V(<p) driving inflation unspecified except to assume 
that it is sufficiently flat in the usual sense; that is e <C 1, 
\rj\ -C 1 where 



M 2 



T7 , V 



Ml — 
v y 



(4) 



are the usual slow roll parameters. 

2 

The coupling ^-(0 — </>o) 2 X 2 in © is introduced to 
ensure that the iso-inflaton field can become instanta- 
neously massless at some point 4> = 4>o along the inflaton 
trajectory (which we assume occurs during the observ- 
able range of e-foldings of inflation) . At this moment x 
particles will be produced by quantum effects. (In section 
IIVI we will discuss how particle production and rescatter- 
ing are modified by the inclusion of a mass term [i 2 x 2 for 
the iso-inflaton.) 



where we have defined the characteristic scale 



(7) 



It is straightforward to verify that the simple expression 
([6]) will be a good approximation for (i?|i|) _1 < 0(e,r]) 
which, in most models, will be true for the entire observ- 
able 60 e-foldings of inflation. 

Note that, without needing to specify the background 
inflationary potential V(<j)), we can write the ratio k±/H 
as 



H 
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2ttV 



1/2 



(8) 



1/2 

where =5xl0~ 5 is the usual amplitude of the vac- 
uum fluctuations from inflation. In this work we assume 
fc* > H which is easily satisfied for reasonable values of 
the coupling g 2 > 10~ 7 . In particular, for g 2 ~ 0.1 we 
have k+/H ~ 30. 

The scenario we have in mind is the following. Infla- 
tion starts at some field value <p > 4>o and the inflaton 
rolls toward the point <j> = 4>o- Initially, the iso-inflaton 
field is extremely massive m x 3> H and hence it stays 
pinned in the vacuum, \ = 0, and does not contribute 
to super-horizon curvature fluctuations. Eventually, at 
t = 0, the inflaton rolls through the point <f> — 0o where 
m x = and \ particles are produced. To describe this 
burst of particle production one must solve the following 
equation for the %-particle mode functions in an expand- 
ing universe 



Xk + 3HXk 



k 2 

^ + kit 2 



(9) 



Equations of this type are well-studied in the context of 
preheating after inflation (42| and moduli trapping [33| . 
In the regime fc* > H particle production is fast com- 
pared to the expansion time 1 and one can solve (jTJ]) very 
accurately for the occupation number of the created x 
particles 



(10) 



A. Quantum Production of \ Particles 



Let us first consider the homogeneous dynamics of the 
inflaton field, 4>{t). Near the point (j) — 4>q we can gener- 
ically expand 



(j){t) ^<j>o+vt 



(5) 



1 In the opposite regime, k+ <C H, the field \ wm be light as 
compared to the Hubble scale for a significant portion of inflation. 
In this case it is no longer consistent to treat the background 
dynamics as being effectively single-field, hence the scenario has 
changed considerably. We will not consider this possibility any 
further. 
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Very quickly after the moment t = 0, within a time At ~ 
k^ 1 <C H^ 1 , these produced x particles become non- 
relativistic (m x > H) and their number density starts to 
dilute as a~ 3 . 

Following the initial burst of particle production there 
are two distinct physical effects which take place. First, 
the energetic cost of producing the gas of massive out-of- 
equilibrium \ particles drains energy from the inflaton 
condensate, forcing <f> to drop abruptly. This velocity 
dip is the result of the backreaction of the produced x 
fluctuations on homogeneous condensate <f)(t). The sec- 
ond physical effect is that the produced massive \ parti- 
cles rescatter off the condensate via the diagram Fig. Q] 
and emit bremsstrahlung radiation of light inflaton fluc- 
tuations (particles) . Backreaction and rescattering leave 
distinct imprints in the observable cosmological pertur- 
bations. Let us discuss each separately. 



B. Backreaction Effects 

We first consider the impact of backreaction. This 
effect can be studied analytically using the mean field 
equation 



where the vacuum average is computed following 



(x 2 > = e(t) 



,-3 



(11) 



(12) 



and n v = J ■ d " 



fc 3 is the total number density 



"x ~ J JTP? nk 
of produced x particles. The Heaviside function &(t) in 
(TT2"]) enforces the fact the the backreaction effects become 
important only for t > 0, after the x particles have been 
produced. The factor of a~ 3 in (Ti~2"j) reflects the usual 
volume dilution of non-relativistic matter. 

The solutions of ((TTJ) display the expected behaviour: 
the energetic cost of the production of x particles at t = 
leads to an abrupt dip in the velocity <p 1 momentarily vi- 
olating the smallness of the slow roll parameter <fr/(H<j>). 
Within a few e-foldings of the moments t — 0, the pro- 
duced x particles have become extremely massive and 
have been diluted away by the inflationary expansion of 
the universe. At this time, the inflaton must settle back 
onto the slow roll trajectory, <j> = —V'/(3H). 

Backreaction effects lead to a transient violation of 
slow roll, and hence we expect an associated "ringing" 
pattern (damped oscillations) in the primordial curva- 
ture fluctuations, similar to models with a sharp feature 
in the potential fl^ - fl6l |491-|54T| . This effect can be seen by 
solving the well-known equation for the curvature pertur- 
bation on co-moving hypersurfaces, 1Z, in linear theory: 



n'l + 2—n' k + k 2 n k = o 

z 



(13) 



Here the prime denotes derivatives with respect to con- 
formal time r = J a _1 (i')<ii' and z = a<j>/H. Note that 



(1131) is valid only in the absence of entropy perturbations. 
However, in our case the x field is extremely massive, 
rn 2 , ^> H 2 , for nearly the entire duration of inflation, 
hence the direct iso-curvature contribution to 1Z is negli- 
gible. 

In [l| the coupled system (|11I13[) was solved numeri- 
cally and the expected ringing pattern in the power spec- 
trum Piz(k) — ^p-|7\Lfc| 2 was obtained. (See also Q.) 
This effect is sub-dominant to the rescattering processes 
described in the next subsection, hence we will not pur- 
sue backreaction any further in this work. 



C. Rescattering Effects 

The second physical effect which takes place after the 
quantum production of x particles in the model (JlJ is 
rescattering. This effect was considered for the first time 
in the context of inflationary particle production in [lj. 
Fig. Q] illustrates the dominant process: bremsstrahlung 
emission of long-wavelength 6<f> fluctuations from rescat- 
tering of the produced x particles off the condensate 4>{t). 
The time scale for such processes is set by the microscopic 
scale, k~ l , and is thus very short compared to the ex- 
pansion time, H~ 1 . Moreover, the production of inflaton 
fluctuations 5(p deep in the IR is extremely energetically 
inexpensive, since the inflaton is very nearly massless. 
The combination of the short time scale for rescatter- 
ing and the energetic cheapness of radiating IR 8<f> leads 
to a rapid build-up of power in long wavelength inflaton 
modes: IR cascading. This effect leads to a bump-like 
feature in the power spectrum of inflaton fluctuations, 
very different from the ringing pattern associated with 
backreaction. The bump-like feature from rescattering 
dominates over the ringing pattern from backreaction for 
all values of parameters. 

In [l[ the model © was studied using lattice field 
theory simulations, without neglecting any physical pro- 
cesses (that is to say that full nonlinear structure of the 
theory, including backreaction and rescattering effects, 
was accounted for consistently). However, this same dy- 
namics can be understood analytically by solving the 
equation for the inflaton fluctuations 8<f> in the approxi- 
mation that all interactions are neglected, except for the 
diagram Fig. [TJ The appropriate equation is 



3H5<f> 



< + V,^8<i ) =-g 2 m)-<p ]x 2 (14) 



The solution of (fT4|) may be split into two parts: the 
solution of the homogeneous equation and the particular 
solution which is due to the source term. Schematically 
we have 



6<f>(t, x) = (ty vac (f,x) + 6<f> ICS c(t, x) 

homogeneous particular 



(15) 



The former contribution is the homogeneous solution 
which behaves as 50 V ac ~ H/(2n) on large scales and, 
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physically, corresponds to the usual scale invariant vac- 
uum fluctuations from inflation. The particular solu- 
tion, 8<p Tesc , corresponds physically to inflaton fluctua- 
tions which are generated by rescattering. The abrupt 
growth of x inhomogeneities at t = sources the par- 
ticular solution 8<f) ICSC , leading to the production of in- 
flation fluctuations which subsequently cross the horizon 
and become frozen. 

A detailed analytical theory of equation (fT4|) will be 
the subject of sections IIVI and [V] Here we simply point 
out that the primordial power spectrum in the model ([1]) 
may, to good approximation, be described by a simple 
semi-analytic fitting function Q 



Pn(k) = A s 



A m ( — 



.3/2 



k 

km 



(16) 



where the first term corresponds to the usual vacuum 
fluctuations from inflation (with amplitude A s and spec- 
tral index n s ) while the second term corresponds to the 
bump-like feature from particle production and IR cas- 
cading. The amplitude of this feature (^ir) depends on 
g 2 , while the location (&ir) depends on <fio. 

In [2| the simple fitting function (jTB")) was used to place 
observational constraints on inflationary particle produc- 
tion using a variety of cosmological data sets. Current 
data are consistent with rather large spectral distortions 
of the type ([To]) . Features as large as Ayh/A s ~ 0.1 are 
allowed in the case that &ir falls within the range of scales 
relevant for CMB experiments. A feature of this magni- 
tude corresponds to a realistic coupling g 2 ~ 0.01. Even 
larger values of g 2 are allowed if the feature is localized on 



smaller scales. In [55[ Large Scale Structure forecast con- 
straints were considered for the model (fT6"|) . It was shown 
that, for Air ^ 0.1 Mpc" 1 , the constraint on Air/j4 s will 
be stengthened to the 0.5% level by Planck or 0.1% in- 
cluding also data from a Square Kilometer Array (SKA). 
With a Cosmic Inflation Probe (CIP) similar constraints 
could be achieved for fciR as large as 1 Mpc -1 . 



D. Nongaussianity from Particle Production 

The bump-like feature in P(k), corresponding to the 
second term in (IT!))) , must be associated with a nongaus- 
sian feature in the bispectrum [![. Indeed, it is evident 
already from inspection of equation (TT¥| that the inflaton 
fluctuations generated by rescattering are significantly 
nongaussian; the particular solution of (1141) is bi-linear 
in the gaussian field \- 

Nongaussian statistics have attracted a considerable 
amount of interest recently, owing to their potential as 
a tool for observationally discriminating between the 
plethora of inflationary models in the literature. Al- 
though the simplest single-field slow roll models are 



known to produce negligible (primordial) nongaussian- 
ity [56l - (58| , there are a currently a number of alter- 
native models which may predict an observable signa- 
ture. Examples include models with preheating into 
light fields [HJ, I59W61I . nonlocal inflation I62B . the cur- 
vaton mechanism |63| . multi- field models 64j], construc- 
tions with a small sound speed [65| (such as DBI [66[ 
inflation), trapped inflation Q, the gelaton [67[, models 
with features or rapid oscillations in the inflaton poten- 
tial [U HH, non- vacuum initial conditions [H, l68l - [70l |. 
warm inflation 71], etc. 

Nongaussianity is usually characterized in terms of the 
bispectrum, B^(ki), which is the 3-point correlation func- 
tion of the Fourier transform of the primordial curvature 
fluctuation on uniform density hypersurfaces, £. Explic- 
itly, we define 

(Ck.Ck.Cka) = (2 7 r) 3 ,5( 3 )(k 1 +k 2 + k 3 )B c (fc 4 ) (17) 

where ki = |ki| and (fc is related to the variable IZk ap- 
pearing in (fiU)) as Cfc — —7Zk on large scales k <§; aH . The 
delta function in (|17p reflects translational invariance and 
ensures that B^(ki) depends on three wavenumbers which 
form a triangle: ki + k 2 + k3 = 0. A general bispectrum 
B^(ki) may be characterized by specifying its size (am- 
plitude of B(), shape (whether B^ peaks on squeezed, 
equilateral or flattened triangles) and running (the de- 
pendence of B^ on the size of the triangle) . The various 
nongaussian scenarios discussed above may be classified 
according to the size, shape and running of the bispec- 
trum, see (72j for a more detailed review. 

The nongaussian signature from IR cascading is very 
different from other models, such as the local, equilateral 
or enfolded shapes, which have been studied in the lit- 
erature. IR cascading only influences modes leaving the 
horizon near the moment <fi — cf>o, when particle produc- 
tion occurs, hence we expect the bispectrum to be very 
far from scale invariant (this is also true for the model 
considered in [Hll [52j]). The dominant contribution to 
B^(ki) should peak strongly for triangles with a charac- 
teristic size ~ fciR, corresponding to the location of the 
bump in the power spectrum l|16p . We will estimate the 
shape of the bispectrum from particle production and IR 
cascading in more detail in section IIVI and re- visit this 
issue also in an upcoming publication [73j |. 

The unusual shape and strong scaling properties of the 
bispectrum from particle production makes it difficult to 
compare the magnitude of nongaussianity in this model 
to more familiar bispectra, such as the local shape, which 
are very close to scale invariant. We find it useful to quan- 
tify the magnitude of the nongaussianity in the model (UJ) 
by computing the moments of the Probability Distribu- 
tion Function (PDF), P((), which is the probability that 
the curvature perturbation has a fluctuations of size (. 
These moments carry information about the correlation 
functions of ( integrated over all wavenumbers ki and 
therefore provide a useful tool to compare models with 
very different shape/running properties [7T| . (See also 
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78] for a related discussion and alternative methodol- 
ogy-) 

Let us define the central moments of the PDF as 



(18) 



The n-th cummulant k„ is the connected n-point func- 
tion. For (£) = the first few non- vanishing cummulants 
are: 



K2 
K5 



s-2\ _ 2 
0> 



(C 2 
(C 3 > 

(C 4 )-3(C 2 ) 2 
<C 5 ) - io(C 3 XC 2 ) 



(19) 
(20) 
(21) 
(22) 



It is useful to introduce the dimensionless cummulants, 
defined as 



(23) 



'C 



For a gaussian PDF we have k n = for n > 3, hence 
these quantify departures from gaussian statistics. When 
the nongaussianities are small, |K n >3| <C 1, then the cor- 
rections to -P(C) are well described by the Edgeworth 
expansion: 



p(0 = 



1 



= -C 2 /(2<t?) 



(24) 



where H^{x) = x 3 — 3x is a Hermite polynomial and the 
• • • denotes corrections of order ka, k\ and smaller. Sec 
[77l [79! . [80l | for more details and [8l| for an alternative 
derivation. 



E. Relation to Other Works 

Before proceeding to study the model ^ in detail, it is 
worth commenting on the relationship between our anal- 
ysis and previous works. Trapped inflation 0, |33| is a 
very closely related model that uses multiple bursts of 
particle production, each similar to the event described 
in subsection lll A\ in order to slow the motion of the infla- 
ton on a steep potential. In that case, dissipation which 
results from the backreaction effects discussed in subsec- 
tion III Bl actually dominate over the friction term 3H <fi 
which is due to the inflationary expansion of the universe. 
In contrast, here we assume that the inflaton potential 
V(4>) is sufficiently flat to support slow-roll inflation; see 
equation ((4]). In our scenario, dissipative effects on the 
homogeneous motion of <p(t) due to particle production 
(backreaction effects) are always negligible; see subsec- 
tion III Bl Nevertheless, we expect that many of our re- 
sults and analytical techniques will also be applicable in 
the case of trapped inflation. 

Trapped inflation was not the first model to attempt to 
exploit dissipative effects in order to assist in slowing the 



motion of the inflaton. Another interesting construction 
of this type is warm inflation [24Tj27| . which employs 
a gas of particles in thermal equilibrium. Again, our 
analysis is distinguished from this mechanism since we 
assume that the homogeneous dynamics of the inflaton 
4>(t) are of the usual slow-roll type. 

Our main focus in this work is to determine the ob- 
servational consequences of isolated burst of particle pro- 
duction on cosmological observables such as the spectrum 
and nongaussianity of the primordial fluctluations. In 
this sense, the spirit of our investigation is more similar 
to works such as 043 than to warm inflation or trapped 
inflation. However, unlike those papers, we have con- 
sistently accounted for rescattering effects which provide 
the dominant contribution to observables in the case at 
hand; see subsection lll CI 



III. NONGAUSSIANITY OF THE 
PROBABILITY DISTRIBUTION FUNCTION 

In order to quantify the magnitude of the nongaus- 
sianity generated by particle production, let us now con- 
sider the PDF in the model ([T]). We proceed numeri- 
cally, re-visiting the lattice field theory simulations per- 
formed in [l| using the HLattice code (74|. For illus- 
tration, we assume the standard chaotic inflation model 
V{4>) = to 2 2 /2 with to = 1Q- 6 V8ttM p and (j) = 



3.2v87rAf p . We consider three different choices of cou- 
pling, g 2 — 1, 0.1, 0.01. Our simulations are performed in 
a 512 3 box whose co-moving size is initially ~ 3 times the 
horizon size H . We run our simulations for roughly 3 
e- foldings from the moment when <j) = cf>o, which is more 
than enough to see the feature from IR cascading freeze 
out as an observable, super-horizon density fluctuation. 
Our choice of (f>o ensures that the feature will be frozen-in 
at scales slightly smaller than the current horizon. Note 
that our quantitative results do not depend sensitively 
on the choice of 4>o 1 nor 011 the details of the background 
inflationary potential V (</>); see [l| for further discussion. 

We extract the PDF of 8(f) from our HLattice simu- 
lations by measuring the fraction of the simulation box 
which contains the fluctuation field 8<j) at a particular 
value. Notice that this approach is completely nonpertur- 
bative: it does not rely on the validity of the Edgeworth 
expansion, nor does it assume anything about the size or 
ordering of the cummulants. This procedure implicitly 
puts a IR cut-off at the box size L and a UV cut-off at 
the lattice spacing, A -1 . Since the nongaussian effects 
in our model are strongly localized in Fourier space, our 
quantitative results are largely insensitive to L and A. 

In Fig. [5] we plot our numerical result for the PDF of 
the inflaton fluctuations generated by rescattering and 
IR cascading. In order to make the physics of inflation- 
ary particle production clear, we have subtracted off the 
contribution coming from the usual vacuum fluctuations 
of the inflaton. That is, the PDF in Fig. [5] is associated 
only with the contribution 8(j) vesc in equation (1151) . 
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FIG. 2: The PDF of the inflaton fluctuations generated by 
rescattering and IR cascading, at a series of different values of 
the scale factor, a. The dotted black curve shows a Gaussian 
fit at late times and we have normalized the scale factor so 
that a = 1 at the moment when particle production occurs. 
For illustration, we have chosen g 2 = 0.1 and a standard 
chaotic inflation potential V(4>) = m 2 (f> 2 /2. 



We can understand physically the behaviour of PDF 
plotted in Fig. [2J Shortly after the initial burst of par- 
ticle production the inflaton perturbations 5<fi are ex- 
tremely nongaussian, due to the sudden appearance of 
the source term J oc x 2 in the equation of motion ([141) . 
Very quickly, in less than an e-folding, nonlinear interac- 
tions begin to drive the system towards gaussianity. A 
very similar behaviour has been observed in lattice simu- 
lations of out-of-equilibrium interacting scalar fields dur- 
ing preheating [75|, [7(| . In the case of rescattering during 
preheating, the system will eventually become gaussian 
when the fields thermalize. However, in our case the uni- 
verse is still inflating. As a result, nongaussian inflaton 
fluctuations generated by rescattering are stretched out 
by the quasi-de Sitter expansion and must freeze once 
their wavelength crosses the Hubble scale. Hence, at late 
times the PDF does not become completely gaussian, but 
rather freezes-in with some non-trivial skewness. Within 
a few e-foldings from the moment of particle production 
the time evolution of the PDF has become completely 
negligible. 

In order to characterize the nongaussianity of the ob- 
servable primordial fluctuations, we would like to con- 
struct the PDF for the curvature perturbation (, includ- 
ing both the contributions from the vacuum fluctuations 
of the inflaton and also from rescattering. To this end, 
we construct £ using the naive relation £ = —£5<j) (see 

section [V] for justification) and take into account both 
contributions to 5<f> in equation (1151) . In Fig. [3] we plot 



FIG. 3: The PDF of the total curvature fluctuation, at late 
times (well after all relevant modes have crossed the horizon 
and frozen). The solid black curve is the exact result from 
our HLattice simulations and the dotted red curve is a gaus- 
sian fit. We have also plotted the leading correction to the 
gaussian result in the Edgeworth expansion, given explicitly 
by equation (24). For illustration, we have chosen g 2 = 0.1 
and a standard chaotic inflation potential V((f>) = m 2 (f> 2 /2. 



the full PDF obtained in this manner, evaluated at very 
late times, well after all relevant modes have crossed the 
horizon and become frozen. 

Given our numerical results for the PDF of the to- 
tal observable curvature fluctuation, that is Fig. [3J it 
is straightforward to compute dimensionless cummulants 
([23| such as the skewness (£3), kurtosis (£4) and noltosis 
(£5) for various values of the coupling g 2 . We have sum- 
marized our results in Table HI Note that for g 2 = 0.01 
both K4 and £5 are too small to be measured accurately 
from our HLattice simulations. 

In order to give some sense of the magnitude of the 
nongaussianity from particle production we have also 
computed an "equivalent /jyff 1 " defined by 5k 3 /(18o-£) 
where the variance is = (C 2 ) 1//2 ~ 10~ 9 / 2 . For a 
given g 2 this effective is the magnitude of Jnl 

which would be necessary to reproduce the skewness 
K3 of the IR cascading PDF using a local ansatz C, = 

C 9 + lfNL [C S 2 -(C S 2 >]. 2 

From Table Q] we see that IR cascading during infla- 
tion can generate significant nongaussianity. Even tak- 
ing g 2 = 0.01 (which is compatible with cosmological 



2 Our sign conventions for /jvl are consistent with WMAP |82( . 
See [77j for a discussion of various conventions employed in the 
literature. 
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TABLE I: Moments of the Probability Distribution Function 



g J skewness kurtosis 


5-th moment 


"equivalent" 


(k 3 ) (k 4 ) 


(%) 


rlocal 
JNL 


1 -0.51 0.2 


1.2 


-4500 


0.1 -0.49 -0.1 


1.5 


-4300 


0.01 -0.006 < O(10 -3 ) 


< e>(nn 3 ) 


-53 



data for any choice of 0o 3) we still obtain a skewness 
«3 = —0.006, which is the same value that would be pro- 
duced by a local model with /jvl ~ — 53. This "equiv- 
alent" local nongaussianity is comparable to current ob- 
servational bounds, and is well within the expected accu- 
racy of future missions. This suggests that nongaussian 
features from particle production during inflation might 
be observable for reasonable values of g 2 . 

The "equivalent" /]^£ al values presented in Table U 
must be interpreted with care. We have included this 
information only to give a heuristic sense of the magni- 
tude of nongaussianity in our model. It must be stressed 
that the PDF plotted in Fig. [3] is quite different from 
the analogous result for local-type nongaussianity. For 
example, the value of the kurtosis (and higher moments) 
are different, as is the ordering of the cummulants. More- 
over, we should emphasize that observational bounds on 
InT 1 cannot be directly applied to our model since the 
bispectrum in our case is uncorrelated with the vacuum 
fluctuations and is far from scale invariant. A detailed 
study of the detectability of nongaussianity from particle 
production will be the subject of a upcoming publication 

Depending on the value of <f>o, the model ([l) may lead 
to a variety of observable signatures. As discussed previ- 
ously, </>o controls the location of the feature in the pri- 
mordial power spectrum (fl~6)). Nongaussian effects are 
also localized near the same characteristic scale, fciR. If 
fciK corresponds to scales relevant for CMB experiments, 
then we predict a bump-like feature in the primordial 
power spectrum, P^(k), and an associated feature in the 
bispectrum B^(ki) with an unusual shape (that will be 
discussed in section HV)) . A key question is whether the 
nongaussian feature can be observably large in a regime 
where the power spectrum feature is small enough to be 
compatible with current observations. Preliminary re- 
sults are encouraging: for g 2 = 0.01 the power spec- 
trum is consistent with all observational data while 
the skewness of the PDF is rather large. A detailed in- 
vestigation will require a simple, separable template for 
the bispectrum and will be discussed in a future publica- 
tion [73j . 

On the other hand, we could imagine a scenario where 
the feature from IR cascading shows up on smaller scales, 
relevant for Large Scale Structure (LSS) experiments 
77, 83 85]. In this case our scenario could be probed 
using higher order correlations of LSS probes (such as 
the galaxy bispectrum) or the abundance of collapsed ob- 



jects (or voids). The latter possibility is interesting since 
the cluster/void abundance is determined the tails of the 
PDF and may be insensitive to the detailed shape of the 
bispectrum. Quantitative predictions for observable clus- 
ter/void abudances require the PDF of the evolved den- 
sity field, smoothed on some relevant scale (73|, rather 
than the PDF of the primordial curvature perturbation 
(which is plotted in Fig. [3]). However, we can neverthe- 
less describe the qualitative signatures which should be 
expected. Our model robustly predicts a negative skew- 
ness for both the curvature perturbation, £, and the den- 
sity field, 5p/p. Hence, we should expect a decrease in 
the abundance of the largest collapsed objects and an 
increase in the abundance of the largest voids [13, [H[ . 
Owing to the localized nature of the bispectrum feature, 
we expect that this effect should show up only when the 
density field is smoothed on a scale close to fciR. 

It is worth mentioning that recent weak lensing mea- 
surement of the dark matter mass of the high-redshift 
galaxy cluster XMMUJ2235.3-2557 [H have been con- 
strued as a possible hint of nongaussian initial conditions 
[90j . Unfortunately, our model does not produce the cor- 
rect sign of skewness to explain such observations. 



IV. ANALYTICAL FORMALISM 

In [l| we studied particle production, rescattering and 
IR cascading using nonlinear lattice field theory simula- 
tions. In addition to this numerical approach, a cursory 
analytical formalism was also presented. Here we re- visit 
the analytical analytical theory of particle production, 
rescattering and IR cascading in an expanding universe, 
in order to better understand the results of [1| from a 
physical perspective. 



A. The Prototype Model 

We consider now the theory 
r2 



S 




M 2 | 

-JLR--(d<t>) 2 -V(0) 



H 

2 



{4> - <M V 



(25) 



which differs from our original model ([2]) by the inclu- 
sion of a mass term AC = —p 2 x 2 /2 for the iso-inflaton. 
Such a term is not forbidden by any symmetry and hence 
one typically expects it to be generated by radiative cor- 
rections, even if the iso-inflaton is classically massless at 
(f> = <Pq. The new parameter /i has the effect of reduc- 
ing the efficiency of the particle production effects dis- 
cussed in subsection III Al the time- varying mass of the 
iso-inflaton 
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does not vanish at (f> = <f)Q, but rather reaches a minimum 
value [x 2 , making the adiabaticity condition more difficult 
to violate. A concern is the possibility that radiative 
corrections induce a large /i and suppress the observable 
effects associated with inflationary particle production. 
Indeed, it is well known that fine tuning may be required 
to keep the mass of any scalar field significantly below 
the cut-off scale associated with the validity of the effec- 
tive field theory description (|25|) . Below, we will show 
that that the suppression of x-p&rticle production is not 
significant provided the following condition is satisfied 

« kl (26) 

where k± = ^/g|u|. Depending on how (|25[) is embedded 
within a more complete framework, the constraint (|26[) 
may (or may not) require fine-tuning to satisfy. Below, 
we will show that the condition (|26p is quite naturally 
satisfied for a large number of microscopically realistic 
constructions. 

Our prototype model (|25l) has been choosen to elu- 
cidate the key physics and observational signatures of 
inflationary particle production in a simple framework 
wherein computations are tractable. We expect, how- 
ever, that many of our qualitative results will carry over 
to more complicated scenarios. In particular, one might 
wish to supplement the action (|25[) by its supersymmetric 
(SUSY) completion; see the interesting work [27| for an 
explicit example. Such an embedding has the advantage 
that the flatness of the inflaton potential V(4>) may be 
protected from large radiative corrections coming from 
loops of the x field. Moreover, a SUSY embedding of the 
model (j2"5)) also allows for some control over the quantum 
corrections to the mass scale \i. 

For models obtained from string theory or super- 
gravity (SUGRA), it is natural to have /z of order the 
Hubble scale 3 during inflation 38-40]; hence we expect 
/i 2 ~ H 2 for such models. In that case, the constraint 
([26"1) is automatically satisfied, because k 2 ^> H 2 when- 
ever particle production is fast as compared to the ex- 
pansion time (that is, for reasonable values of the cou- 
pling g 2 > 1CU 7 , which we assume throughout this work). 
Hence, there exists a very large class of realistic micro- 
scopic models in which radiative effects will not spoil the 
observational consequences of inflationary particle pro- 
duction and IR cascading. 

Although the condition for the efficiency of particle 
production - that is equation (|26|) - can be easily satis- 



3 In the context of SUGRA, the finite energy density driving infla- 
tion breaks SUSY and induces soft scalar potentials with curva- 
ture of order V^ ft ~ fi 2 ~ H 2 In the case of string theory, 
many scalars are conformally coupled to gravity [4QJ through an 
interaction of the form SC = — j^Rx 2 where the Ricci scalar 
is R ~ H 2 during inflation. More generally, any non-minimal 
coupling SC = —^Rx 2 between gravity and the iso-infiaton will 
induce a contribution of order H to the effective mass of x, as 
long as £ = 0(1). 



lied for models coming from string theory or SUSY, we 
prefer to remain agnostic regarding how the prototype ac- 
tion ()25[) is embedded within a more complete framework. 
Throughout our analysis we will keep the inflaton poten- 
tial V(cj>) and the iso-inflaton mass parameter \x more-or- 
less arbitrary (we assume that the slow roll conditions are 
satisfied, and also that /i 2 > 0). This phenomenological 
approach is not different from the philosophy that is em- 
ployed in the majority of work on inflationary cosmology, 
since the slow roll conditions (Q| may be sensitive to UV 
physics whose detailed form is often not specified. The 
question of how the model (|23|) . with a given choice of 
V(<j>) and /x, arises from some complete model of particle 
physics is interesting. However, this question it is not 
the main focus of the current investigation. We refer the 
reader to 2] for several example microscopic embeddings 
within string theory and also SUSY (see also Q). 

Let us now proceed to develop an analytical formalism 
to study inflationary particle production in the model 
(|25|) . The equations of motion that we wish to solve are 

-□^ + U'(^+. 9 2 (0-0o)x 2 = (27) 

-□x+ [n 2 + g 2 {4>-M 2 ]x = o (28) 

where □ = s^VV" is the covariant d'Alembertian. It 
will be useful to work with conformal time r, related to 
cosmic time t via adr = dt. In terms of conformal time 
the metric takes the form 

ds 2 = -dt 2 + a 2 (t)dx- dx 

= a 2 (r) [-dr 2 + dx • dx] (29) 

We denote derivatives with respect to cosmic time as / = 
dtf and with respect to conformal time as /' = d T f. The 
Hubble parameter H = a/a has conformal time analogue 
% = a! J a. For an inflationary (quasi-de Sitter) phase 
{H = const) one has 



to leading order in the slow roll parameter eCL 

As discussed in section HH the motion of the homoge- 
neous inflaton <j>(t) leads to the production of a gas of 
X particles at the moment t — when <p — <j>Q. The 
first step in our analytical computation is to describe 
this burst of particle production in an expanding uni- 
verse. Following the initial burst, both backreaction and 
rescattering effects take place. Our formalism will focus 
on the latter effect, which has been shown to be much 
more important [ij. 



B. Particle Production in an Expanding Universe 

The first step in our scenario is the quantum mechan- 
ical production of x-particles due to the motion of <fi. To 
understand this effect we must solve the equation for the 
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X fluctuations in the rolling inflaton background, 
proximating <f> = 4> + vt equation (|28[) gives 



V 2 



+ [if + kit 2 ] X = 



Ap- 



(31) 



where fc* = y5|u]. We remind the reader that k* H 
for reasonable values of the coupling, see equation ((5J. 

The flat space analogue of equation (|3ip is very well 
understood from studies of broad band parametric res- 
onance during preheating [42| and also moduli trapping 
at enhanced symmetry points (33j . One does not expect 
this treatment to differ significantly in our case since both 
the time scale for particle production At and the char- 
acteristic wavelength of the produced fluctuations A are 
small compared to the Hubble scale. Hence, we expect 
that the occupation number of produced x particles will 
not differ significantly from the flat-space result, at least 
on scales k > H . Furthermore, notice that the x is 
extremely massive for most of inflation. Indeed, even in 
the case /i 2 = we have 



"*x 
H 2 



kte 

H 2 



Since fc* 3> H it follows that m 2 



(32) 



3> H 2 , except in a 



tiny interval if|Af| ~ (H/k+) 2 which amounts to roughly 
10~ 3 e- foldings for g 2 ~ 0.1. Therefore, we do not ex- 
pect any significant fluctuations of x to be produced on 
super-horizon scales k < H. (Allowing for /i 2 ^ only 
strengthens this conclusion.) 

Let us now consider the solutions of equation (|3"Tj) . We 
work with conformal time r and write the Fourier trans- 
form of the quantum field x as 



x{t, x) 



d 3 k 



(2tt) 3 / 2 o(r) 



(33) 



Note the explicit factor of a" 1 in (|3"3")l which is introduced 
to give a canonical kinetic term. The q-number valued 
Fourier transform £k( T ) c & n be written as 



£k( r ) =OkX*(T) + al k x*(r) 



(34) 



where the annihilation/creation operators satisfy the 
usual commutation relation 



a k ,4< =<5 (3) (k-k') 



(35) 



and the c-number valued mode functions Xk obey the 
following oscillator-like equation 



x'l(r)+^ k (r)Xk(r)=0 
The time-dependent frequency is 



(36) 



k 2 + a 2 m 2 (r) 



a 
a 



k z + — 



(37) 



where 



2 / \ 2 , 2 

m (r) = A* +.9 



= ^ + ^V(r) 



(38) 



is the time-dependent effective mass of the x particles 
and 



(39) 



is the usual cosmic time variable. We have arbitrarily set 
the origin of conformal time so that r = —1/H corre- 
sponds to the moment when <f> — (f>Q. 

It is useful to define the occupation number riu of the 
X of particles with momentum k, defined as the energy 
of the mode ||x'fc| 2 + 5 w felXfc| 2 divided by the energy w fe 
of each particle. Explicitly, we define 



2 



\X k 



l 12 



\Xk\ 



(40) 



where the term — ^ comes from extracting the zero-point 
energy of the linear harmonic oscillator (see [42| for a 
review). Our definition (|40)) coinicides with the usual 
notion of particle number in the asymptotic adiabatic 
regimes (|i| >, A:" 1 ). During the very brief non-adiabatic 
period (|t| <C k^ 1 ) our result coincides with the usual no- 
tion of quasi-particle number, obtained by instantaneous 
diagonalization of the Hamiltonian 

Let us now try to understand analytically the be- 
haviour of the solutions of (f36|) . At early times t <C — k~ 1 , 
the frequency ujk varies adiabatically 



a) 7 



< 1 



(41) 



In this in-going adiabatic regime the modes Xfc are not 
excited and the solution of (fBTij) are well described by the 
adiabatic solution Xfc( T ) = /fc( T ) where 



fk(r) 



1 



exp 



i / dr 'wfe(r) 



(42) 



We have normalized (|42[) to be pure positive frequency 
so that the state of the iso-inflaton field at early times 
corresponds to the adiabatic vacuum with no x particles. 
(Inserting (|4*2"j) into P0|) one finds n k = for the adiabatic 
solution, as expected.) 

The adiabatic solution ([42"T) ceases to be a good approx- 
imation very close to the moment when <p — <po, that 
is at times \t\ < k^ 1 . In this regime the adiabaticity 
condition (|4ip is violated for modes with wave-number 
H £ k <, \Jk 2 — [l 2 and x particles within this momen- 
tum band are produced. During the non-adiabatic regime 
we can still represent the solutions of in terms of the 
functions /fc(r) as 



Xk(r) = a k (r)f k (r)+I3 k (r)mr) 



(43) 
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This expression affords a solution of (1361) provided the 
time-dependent Bogoliubov coefficients obey the follow- 
ing set of coupled equations 



&(t) = ^\exp 
2uj k (T) 



+2i J dr'uj k (T') 
-2i f dT'oj k (r') 



Pk(r) (44) 
a k (r) (45) 



The Bogoliubov coefficients are normalized as |afc| 2 — 
\Pk\ 2 = 1 and the assumption that no x particles are 
present in the asymptotic past 4 fixes the initial conditions 
a k = 1, Pk — f° r t ~ * —oo. This is known as the 
adiabatic initial condition. 

From the structure of equations (|44I45|) it is clear that 
violations of the condition (|4ip near t = leads to a 
rapid growth in the \/3k\ coefficient. The time variation 
of 0k can be interpreted as a corresponding growth in the 
occupation number. Inserting (|43|) into (|40|) we find 



n k 



\Pk 



(46) 



At late times (t > fc^ 1 ) adiabaticity is restored and 
the growth of n k = \Pk\ 2 must saturate. By inspection 
of equations (|44l45p we can see that the Bogoliubov co- 
efficients must tend to constant values in the out-going 
adiabatic regime. Therefore, within less than an e-folding 
from the moment of particle production the solution Xk 
of equation (|36[) can be represented as a simple super- 
position of positive frequency fk modes and negative fre- 
quency /£ modes. Our goal now is to derive an analytical 
expression for the modes Xk which is valid in this out- 
going adiabatic region. 

First, we seek an expression for the Bogoliubov coeffi- 
cients a k , fik in the out- going adiabatic regime t £ fc^ 1 . 
From (|44I45I) it is clear that the value of the Bogoliubov 
coefficients at late times can depends only on dynamics 
during the interval \t\ < fc" 1 where the adiabaticity con- 
dition (14T1) is violated. This interval is tiny compared to 
the expansion time and we are justified in treating a(r) 
as roughly constant during this phase. Hence, it follows 
that the flat space computation of the Bogoliubov coeffi- 
cients [H, must apply, at least for scales k > H . To 
a very good approximation we therefore have the well- 
known result 



a k = VT+ e -nt* 2 /k 2 te -irk 2 /k 

Pk = —ie 



7rp 2 /(2feJ) e -7rfe 2 /(2feJ) 



(47) 
(48) 



in the out-going adiabatic regime. Equation (I48[) gives 
the usual expression 5 for the co-moving occupation num- 
ber of particles produced by a single burst of broad-band 



parametric resonance: 



(49) 



Comparing equations (|49|) and ([T0|) we see that the mass 
parameter /i for the iso-inflaton has the effect of sup- 
pressing the number density of produced x particles by 
an amount e " M ' * < 1. This suppression reflects the re- 
duced phase space of produced particles: the adiabaticity 
condition is violated only for modes with k < y/k^ — jjp. 
Notice that the suppression of x particle production is 
negligible when [i 2 <C fc 2 , precisely the condition (j2"6")l 
that was alluded to earlier. For the remainder of this 
work we will assume that fi fc*, since in the opposite 
regime the observational signatures of particle produc- 
tion effects are exponentially suppressed. 

Next, we seek an expression for the adiabatic solution 
/fe(r) in the out-going regine t £ fc" 1 . We assume fi £ fc* 
and also fo cus on th e interesting region of phase space, 
H < k < \Jk1 — /i 2 . In this case, the adiabatic solution 
(1421) is very well approximated by 



fk(r) 



1 



\klt\r) 



(50) 



where t{r) is defined by ([59^1 . It is interesting to note 
that equation (|50|) is identical to the analogous flat-space 
result [l(, except for the factor of a -1 / 2 . Taking into 
account also the explicit factor of a -1 in our definition of 
the Fourier transform (|3"3")1 we recover the expected large- 
scale behaviour for a massive field in de Sitter space, that 
is x ~ a 3 ^ 2 - This dependence on the scale factor is easy 
to understand physically, it simply reflects the volume 
dilution of non-relativistic particles: p x ~ m 2 % 2 ~ a -3 . 
Notice that the parameter /x does not appear in (|50l) . 
This is so because, for the time-varying mass of the x 
field, equation (|3"5|) . is dominated by the interaction term 
when k+\t\ 1 and /i fc*. 

Finally, we arrive at an expression for the out-going 
adiabatic x modes which is accurate for interesting scales 



H <k < y/kY 



Putting together the results ([50)1 



and (l43l) along with the well-known expressions (|4TI48|) 
we arrive at 



Xk 



( T ) S Vl + e-^ 2 /fe? e --fe 2 /fc?_ 

n 1 



1 



(51) 



valid for t > k^ 1 . Equation (|5T1) is the main result of 
this subsection. We will now justify that this expression 
is quite sufficient for our purposes. 



4 This assumption is justified since any initial excitation of \ would 
have been damped out exponentially fast by the expansion of the 
universe. 

5 Our result for the Bogoliubov coefficients is consistent with |33H . 



In that work /i was interpreted as an "impact parameter" for 
the motion of the modulus, whereas in our work we interpret 
this as a "bare" mass term. This distinction has no impact on 
the result for the occupation number because, in both cases, the 
parameter appears in the same way in the equation of motion 
for the fluctuations of \- 
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For modes deep in the UV, k > k±, our expression (I5T1) . 
is not accurate. 6 However, such high momentum parti- 
cles are not produced, the condition (j4Tj) is always satis- 
fied for k k+. Note that the absence of particle pro- 
duction deep in the UV is built into our expression (|5ip: 
as k — > oo this function tends to the vacuum solution 
Xk -> Ik- 

Our expression (|51l) is also not valid deep in the IR, 
for modes k < H. To justify this neglect requires some- 
what more care. Notice that, even very far from the 
point tfi = cf>o long wavelength modes k ^ H should 
not be thought of as particle-like. The large-scale mode 
functions are not oscillatory but rather damp exponen- 



tially fast as x 



~ 3 / 2 . Hence, even if we started with 



some super-horizon fluctuations of x at the- beginning of 
inflation, these would be suppressed by an exponentially 
small factor before the time when particle production oc- 
curs. Any super-horizon fluctuation generated near t = 
would need to be exponentially huge to overcome this 
damping. However, resonant particle production during 
inflation does not lead to exponential growth of mode 
functions. 7 

To verify explicitly that there is no significant effect for 
super-horizon fluctuations let us consider solving equa- 
tion (f5T|) neglecting gradient terms. The equation we 
wish to solve, then, is 



kH 1 



9 , 
4 H 



(a 3 / 2 X ) = (52) 



(For simplicity we take [i — and e — for this para- 
graph, however, this has no effect on our results.) The 
solution of this equation may be written in terms of 
parabolic cylinder functions D v (z) as 



X{t, x) 



1 



7,3/2 



[(i + *)M 



+ C 2 D_ 19H 2. [(-l+t)M] (53) 



For our purposes the precise values of the coefficients C\ , 
Ci are not important. Rather, it suffices to note that for 
k*\t\ >, 1 the function (|53|) behaves as 



X(*,x) - \t\- 1/2 e~ 3Ht/2 x [oscillatory] 



(54) 



6 The expression (1 50 1) for the adiabatic modes fk is not valid at 
high momenta where u>^ = k. 

7 In this regard our scenario is very different from preheating at 
the end of inflation. In the latter case the inflaton passes many 
times through the massless point m x = and there are, cor- 
respondingly, many bursts of particle production. After many 
oscillations of the inflaton field, the \ particle occupation num- 
bers build up to become exponentially large and, averaged over 
many oscillations of the background, the \ mode functions grow 
exponentially. However, in our case case there is only a single 
burst of particle production at t = 0. The resulting occupation 
number ||10|I is always less than unity and the solutions of l|36|l 
never display exponential growth. 



This explicit large-scale asymptotics confirms our previ- 
ous claims that the super-horizon fluctuations of x damp 
to zero exponentially fast, as a -3 / 2 ~ e ~3Ht/2^ ^ g 
cussed previously, this damping is easy to understand 
in terms of the volume dilution of non-relativistic parti- 
cles. (If we had included the parameter /^ 2 in equation 
(|52|) then our conclusions would only be strengthened, 
since this parameter has the effect of making the iso- 
inflaton even more massive.) We can also understand the 
power-law damping that appears in (|54|) from a phys- 
ical perspective. The properly normalized modes be- 
have as a?l 2 x ~ ^k 1 ^ 2 wm lc on large scales we have 
lOk ~ \ m x\ ~ \A- Hence, the late-time damping factor 
i^ 1 / 2 which appears in (|54|) reflects the fact that the x 
particles become ever more massive as <f> rolls away from 
the point 4>q. 

Finally, it is straightforward to see that the function 
(fSU)) does not display any exponential growth near t = 0. 
Hence, we conclude that there is no significant genera- 
tion of super-horizon x fluctuations due to particle pro- 
duction. 8 

In this subsection we have seen that the quantum pro- 
duction of x particles in an expanding universe proceeds 
very much as it does in flat space. This is reasonable since 
particle production occurs on a time scale short compared 
to the expansion time and involves modes which are in- 
side the horizon at the time of production. 



Inflaton Fluctuations 



In section TlV Bl we studied the quantum production of 
X particles which occurs when (f> rolls past the massless 
point 4> = 4>q. Subsequently, there are two distinct phys- 
ical processes which take place: backreaction and rescat- 
tering. As we discussed, the former effect has a negligible 
impact of the observable spectrum of cosmological per- 
turbations and may be neglected. 

In this subsection we study the rescattering of pro- 
duced x particle off the inflaton condensate. The dom- 
inant process to consider is the diagram illustrated in 
Fig. [U corresponding to bremsstrahlung emission of 5<p 
fluctuations (particles) in the background of the external 
field. (There is also a sub-dominant process of the type 
XX ~^ SipScj) which is phase space suppressed.) Taking 
into account only the rescattering diagram illustrated in 
Fig. [1] is equivalent to solving the following equation for 



8 This is strictly true only in the linearized theory. It is possible 
that x particles are generated by nonlinear effects such as rescat- 
tering. However, even such second order \ fluctuations will be 
extremely massive compared to the Hubble scale and must there- 
fore suffer exponential damping a -3 / 2 on large scales. 
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the q-number inflaton fluctuation 



V 2 

md t - — 

a 2 - 



= -g 2 m)-Mx 2 (55) 



where we have introduced the notation : 



V M for the 



inflaton effective mass. (Note that we are not assuming 
a background potential of the form m 2 <^ 2 /2, only that 
V<t>4> ^ in the vicinity of the point </> = </>o-) 

Equation ([55]) may be derived by noting that ([2]) gives 
an interaction of the form g 2 (ip—(f>Q)S(f>x 2 between the in- 
flaton and iso-inflaton, in the background of the external 
field 4>(t). Equivalently, one may construct this equation 
by a straightforward iterative solution of (|2"7)l . 

We work in conformal time and define the q-number 
Fourier transform £ k (r) of the inflaton fluctuation anal- 
ogously to (|3"3"|) : 



d 3 k 



(2tt) 3 / 2 a(r) 



(56) 



(To avoid potential confusion we again draw the attention 
of the reader to the explicit factor a -1 in our convention 
for the Fourier transform.) The equation of motion (|55[) 
now takes the form 



d 2 + k 2 + a 2 m 2 



-gkla{r)t{T) 



a' 
a 
d 3 k 



(2tt) 3 / 



a&£-i» (57) 



The solution of ([ST)) consists of two parts: the solution of 
the homogeneous equation and the particular solution 
which is due to the source. The former corresponds, 
physically, to the usual vacuum fluctuations from infla- 
tion. On the other hand, the particular solution corre- 
sponds physically to the secondary inflaton modes which 
are generated by rescattering. 



D. Homogeneous Solution and Green Function 



We consider first the homogeneous solution of (|57[). 
Since the homogeneous solution is a gaussian field, we 
may expand the q-number Fourier transform in terms 
of annihilation/creation operators 6k, b^ and c- number 
mode functions <j>k{j) as 



^(r) = &k^(r) + & t _ k ^(r) 



(58) 



Here the inflaton annihilation/creation operators b^, 6 k 
obey 



b k ,bl 



= *< 3 >(k-k') 



(59) 



and commute with the annihilation/ creation operators of 
the x-ficld: 



[ok, h 



ak,6 k / =0 



(60) 



Using ([30| and ((4]) it is straightforward to see that the 
homogeneous inflaton mode functions obey the following 
equation 



5 2 </> fe + 



= 



where we have defined 



(61) 



(62) 



The properly normalized mode function solutions are well 
known and may be written in terms of the Hankel func- 
tion of the first kind as 



</>k(r) 



-<r#«(-fcr) 



(63) 



This solution corresponds to the usual quantum vacuum 
fluctuations of the inflaton field during inflation. 

In passing, let us compute the power spectrum of the 
quantum vacuum fluctuations from inflation. Using the 
solutions (l63l) we have 



(64) 



p; ac (fc) 



fc 3 




2 H 2 




2tt 2 


a 


~ (2^) 2 





on large scales k <C aH. The explicit factor of a~ 2 in 
(|64l) appears to cancel the a -1 in our definition of the 
Fourier transform (|56p . The spectral index is 



n s - 1 = 3 - 2v = 2?/ - 2e 



(65) 



using (PI) . 

Given the solution ([63t of the homogeneous equation, 
it is now trivial to construct the retarded Green function 
for equation (|57[) . This may be written in terms of the 
free theory mode functions ([6"3"]) as 

G fc (T-T') - i6(r-r')[ ^(r)^(r') - ^(r)^(r') ] 



ITT 



9(r 



H^(-kr)HP(-kr')* 



HP(-kr)*HP(-k^) 



(66) 



E. Particular Solution: Rescattering Effects 



We now consider the particular solution of (|57[). This 
is readily constructed using the Green function (|66[) as 



9% 



(2tt) 3 / 2 



(67) 

dr'd 3 k'G k (r - r') a(r')t(r') (r') 



Notice that this particular solution is statistically in- 
dependent of the homogeneous solution (j58[) . In other 
words, the particular solution can be expanded in terms 
of the annihilation/creation operators ak,a k associated 
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with the x held, whereas the homogeneous solution is 
written in terms of the annihilation/creation operators 
6k, associated with the infiaton vacuum fluctuations. 
These two sets of operators commute with one another. 

We will ultimately be interested in computing the n- 
point correlation functions of the particular solution (|68| . 
For example, carefully carrying out the Wick contrac- 
tions, the connected contribution to the 2-point function 
is 

(^^W) = ^ (3) (k 1+ k 2 ) 

x / dT l dT"a{T l )a{T")t{T l )t{r")G kl (r - r')G fea (r - r") 



As a step towards renormalizing the 2-point correlation 
function of infiaton fluctuations from rescattering (|68p . 
let us first consider the simpler problem of renormalizing 
the 2-point function of the gaussian field X . Wc defined 
the renormalized 2-point function in momentum space as 
follows: 



(4 x 1 (*i)4 x 2 (i 2 )> r 



-<&(*i)&(*a)>i 



(71) 



In ([7T|) the quantity (ti)£% (ta))in is the contribution 
which would be present even in the absence of particle 
production, computed by simply taking the solution (|43[) 
with ak = 1, 0k = 0. Explicitly, we have 

x / d 3 k> Xkl ^(T') x t 1 _ k ,(T") Xk ,(T>) x :,(T") (68) (til(h)$ 2 (h)) in = ^(kx +k 2 )f kl (t 1 )t k2 (t 2 ) (72) 



The power spectrum of 5(f> fluctuations generated by 
rescattering is then defined in terms of the 2-point func- 
tion in the usual manner 



2tt 2 



<^(*)^(r)>=«5 (3) (k + k') ir a(r) 2 P, 



;>rcsc 

4> 



(69) 



(The explicit factor of a 2 in the definition (|69|) appears 
to cancel the factor of a -1 in our convention for Fourier 
transforms ([56]) .) 

The total power spectrum is simply the sum of the 
contribution from the vacuum fluctuations (|64j) and the 
contribution from rescattering ()69[) : 



p (fc) = p; ac (fc) + p^ sc (fc) 



(70) 



There are no cross-terms, owing to the fact Ok and 
commute. 



F. Renormalization 

We now wish to evaluate the 2-point correlator (|68[) . In 
principle, this is straightforward: first substitute the re- 
sult (|51[) for the Xk modes and the result (|66|) for the 
Green function into (|68|) . next evaluate the integrals. 
However, there is a subtlety. The resulting power spec- 
trum is formally infinite. Moreover, the 2-point corre- 
lation function (|68p receives contributions from two dis- 
tinct effects. There is a contribution from particle pro- 
duction, which we are interested in. However, there is 
also a contribution coming from quantum vacuum fluc- 
tuations of the X field interacting non-linearly with the 
infiaton. The latter contribution would be present even in 
the absence of particle production, when a k = 1, Pk — 0. 

In order to isolate the effects of particle production on 
the infiaton fluctuations, we would like to subtract off 
the contribution to the 2-point correlation function 
which is coming from the quantum vacuum fluctuations 
of X . This subtraction also has the effect of rendering the 
power spectrum (|69[) finite, since it extracts the usual UV 
divergent contribution associated with the Minkowski- 
space vacuum fluctuations. 



where f k are the adiabatic solutions (|42l) . 

To see the impact of this subtraction, let us consider 
the renormalized variance for the iso-inflaton field, ( X 2 )- 
Employing the prescription (|71[) we have 



(X 2 (r,x)) re 



d 3 k 



(27r) 3 a 2 (r 



\Xk{rt 



2Wfc(T) 



M 



(73) 



where 8m is the contribution from the Coleman- Weinberg 
potential. This proves that our prescription reproduces 
the scheme advocated in [33] . The renormalized variance 
(1731 is finite and may be computed explicitly using our 
solutions (j5"T) . We find 



where 



(X 2 (*,x)) r 



d 3 k 



n x a 



9- 



(2tt>- 



(74) 



(75) 



is the total co-moving number density of produced X par- 
ticles. The result (|74|) was employed in [7( to quantify the 
effect of backreaction on the infiaton condensate in the 
mean field treatment (|11[) . Hence, the renormalization 
scheme (1711) was implicit in that calculation also. 

At the level of the 2-point function, our renormaliza- 
tion scheme is tantamount to assuming that Coleman- 
Weinberg corrections are already absorbed into the def- 
inition of the infiaton potential, V{<f>). In general, such 
corrections might steepen V{<f>) and spoil slow-roll infla- 
tion. Here, we assume that this problem has already been 
dealt with, either by fine-tuning the bare infiaton poten- 
tial or else by including extended SUSY (which can min- 
imize dangerous corrections). See also [33| for a related 
discussion. Note, also, that our renormalization proce- 
dure is equivalent to the quasi-particle normal ordering 
scheme described in (9l| . 

Having established a scheme for remormalizing the 2- 
point function of the gaussian field X , it is now straight- 
forward to consider higher order correlation functions. 
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We simply re-write the 4-point function as a product of 
2-point functions using Wick's theorem. Next, each Wick 
contraction is renormalized as ([71]) . Applying this pre- 
scription to (|68p amounts to 



2_g^ 

(2tt) 3 



^(k!+k 2 ) 



x / dr / dr / 't(Ty(T / ')a(T / )a(T / ')G kl (r - T r )G k2 (r - t' 1 ) 



x\Xk>(T')xi>(T")-f«(t')f* k ,(T")] (76) 
where fk(t) are the adiabatic solutions defined in (|42"j). 
I 



G. Power Spectrum 



We are now in a position to compute the renormalized 
power spectrum of inflation fluctuations generated by 
rescattering, P^, csc (k). We renormalize the 2-point cor- 
relator of the inflaton fluctuations generated by rescatter 
according to (|76j) and extract the power spectrum by 
comparison to (|69|) . We have relegated the technical de- 
tails to Appendix A and here we simply state the final 
result 



P^ c (k) 



16tt 5 



2V2 



(/ 2 (fc,T) 2 + |/i(fc,r)| 2 ) 



2V2 



{l 2 (k,T) 2 -ne[h(k,T)]) 



8 ^ c -^ l j 2 /(2kl) e ^k 2 /(Zkl) , 4^ -57rM 2 /(2^)„-37r 



3^ 



( „..p 2 /(2fe 2 ) e -37rfc 2 /(5fe 2 ) 



Im [/ 1 (fc,r)7 2 (fc,r)] 



(77) 



where the functions 1%, I2 are the curved space general- 
ization of the characteristic integrals defined in [1| . Ex- 
plicitly we have 

h(k,T) = t-t f dT'G k (T- T 'y k * t2 ^ (78) 

h(k,T) = -L / dr'G^T-T') (79) 

The characteristic integral I2 can be evaluated analyti- 
cally, however, the resulting expression is not particularly 
enlightening. Evaluation of the integral I\ requires nu- 
merical methods. More details in Appendix A. Equation 
([77)) is the main result of this section. 

To test our analytical formalism, let us compare the 
result (|77|) with the out-put of fully nonlinear HLattice 
simulations. In Fig.2]we plot our results for P™ sc (k) as a 
function of k, for several time steps in the evolution. We 
have normalized P™ sc (k) to the amplitude of the usual 
vacuum fluctuations from inflation, P™ c (fc) - H 2 /(2ir) 2 . 
This figure illustrates the final stages of IR cascading; we 
see the peak of the bump-like feature slide to k ~ e _3 £;*, 
at which point the associated mode functions S(f>k have 
crossed the horizon and become frozen. At later times in 
the evolution the peak of the feature and also the IR tail 
(~ k 3 ) remain fixed. Modes associated with the UV end 
of the spectrum are still inside the horizon and continue 
to evolve as 5(fik ~ a_1 j which explains the damping of 



the k > e~ 2 k+ part of the spectrum. At late times, the 
shape of the feature that is frozen outside the horizon can 
be very well approximated by the semi-analytic fitting 
function (fTo]l . 

The agreement between our analytical formalism and 
the exact numerical results is quite evident from Fig. @] 
and provides a highly nontrivial check on our calculation. 



H. The Bispectrum 

So far, we have shown how to compute analytically the 
power spectrum generated by particle production, rescat- 
tering and IR cascading in the model @. We found that 
IR cascading leads to a bump-like contribution to the 
primordial power spectrum of the inflaton fluctuations. 
However, this same dynamics must also have a nontrivial 
impact on nongaussian statistics, such as the bispectrum. 
Indeed, it is already evident from our previous analysis 
that the inflaton fluctuations generated by rescattering 
may be significantly nongaussian. From the expression 
([55]) we see that the particular solution (due to rescat- 
tering) is bi-linear is the gaussian field x- 

We define the bispectrum of the inflaton field fluctu- 
ations in terms of the three point correlation function 
as 

^UUi(r))^(2n) 3 a 3 (r)S(k 1+ k 2 +k 3 )B^h) (80) 
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FIG. 4: The power spectrum of infiaton modes induced by 
rescattering. (normalized to the usual vacuum fluctuations) 
as a function of ln(fe/fc*), plotted for three representative time 
steps in the late-time evolution. For each time step we plot 
the analytical result (the solid line) and the data points ob- 
tained using lattice field theory simulations (diamonds). The 
agreement between these two independent results is evident. 
For illustration, we have set /i 2 = 0. 



The factor o? appears in (|80|) to cancel the explicit factors 
of a -1 in our convention (|56[) for the Fourier transform. 
It is well-known that the nongaussianity associated with 
the usual quantum vacuum fluctuations of the infiaton is 
negligible [56j-|58|, therefore, when evaluating the bispec- 
trum (|5D|) we consider only the particular solution (|B"5)) 
which is due to rescattering. Carefully carrying out the 
Wick contractions, we find the following result for the 
renormalized 3-point function 



(t£AA a (T)) Ien = -g^ a(k x + k 2 + k 3 ) I] / dntirMr^ir n) 



X / d 3 p [Xk 1 -p(Tl)Xkl-p('^) ~ /fcl-p(n)/fc 1 _p(T2)] [Xk 3 +p(T 2 )x*k 3 +p(T3) ~ fk :i +p(T2)fk 3+p (T 3 )] 

x [x P (ri)Xp(r 3 ) - /p(ti)/*(t 3 )] 
+ (fc 2 O k 3 ) 



(81) 



where the modes Xfc are defined by (j43|) and fk are the 
adiabatic solutions (l4"2l . On the last line of (I5T1) we have 
labeled schematically terms which are identical to the 
preceding three lines, only with fc 2 and k% interchanged. 
One may verify that this expression is symmetric under 
interchange of the momenta ki by changing dummy vari- 
ables of integration. 



I. Estimating the Shape of the Bispectrum 

It is straightforward (but tedious) to plug the expres- 
sions (f50|) and (f51"T) into (|81"|) and evaluate the integrals. 
The resulting expression is extremely cumbersome and 
not particularly enlightening. We are interested here in 
extracting some information about the shape of the bis- 
pectrum BAki). For this purpose, it suffices to work in 
the flat-space limit, H 0. This will give a reasonable 
qualitative picture of the full result since the entire pro- 



cess of IR cascading occurs over a time scale somewhat 
shorter than the expansion time. In [l[ this same ap- 
proximation was employed to study the power spectrum 
from IR cascading and was found to reproduce the H ^ 
results to good accuracy. 

A detailed calculation of S^(fcj) has been relegated to 
appendix B. Here we simply provide a representative con- 
tribution, in order to give a rough sense of the qualitative 
behaviour: 



Tfc?/(3feJ) 



1 



+ m 2 t 



(82) 

for some constant C. This expression captures some of 
the qualitative features of the full result, in particular 
the dynamical cascading of nongaussianity into the IR 
to generate a localized bispectrum feature. It should be 
stressed that (1821) is a heuristic estimate and not a fit- 
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ting function nor a systematic approximation to the full 
result. Hence, equation (|52l should not be used to make 
quantitative predictions of any kind. 

As anticipated, our expression for B$(ki) peaks only 
over when all wavenumbers are close to the characteristic 
scale corresponding to the location of the bump in (fT6| . 
Therefore, particle production and IR cascading leads to 
a localized nongaussian feature in the bispectrum, rather 
than the nearly scale-invariant signatures that are usually 
considered. We will discuss the phenomenology of this 
new type of nongaussianity in a forthcoming publication 
& 

Now we would like to attempt to characterize the shape 
of the nongaussianity from particle production and IR 
cascading. To this end we define a "shape function" S(ki) 
as follows 



S(ki) = N- 1 (k 1 k 2 k 3 ) 2 B 4 ,(k l 



(83) 



where N is a normalization factor which will not con- 
cern us. 9 The function S(ki) has the advantage that the 
strong fc 6 running of the bispectrum is extracted. Hence, 
any residual scaling behaviour displayed by S(kj) must 
be a result of nonlinear interactions; see also [5ll. |52||. 

Symmetry of the bispectrum under permutations of 
momenta implies that we can focus only on the region 
k\ > k 2 > k 3l to avoid counting the same configuration 
twice. Moreover, the triangle inequality implies that 1 — 
jr- < Therefore we can completely specify the shape 
of the bispectrum for a given size of triangle k by plotting 
S(k,kx2,kx3) in the region X3 < x 2 < 1 and 1 — x 2 < 
X3. (See also [93j.) Because our bispectrum is very far 
from scale-invariant, it follows that this shape function is 
sensitive to the choice of k. Therefore, in Fig.[5]we choose 
several representative choices: ln(fc/fcbum P ) = —1, 0, 1, 2. 

We see that a rich array of shape are possible: for 
k fcbump the bispectrum is qualitatively similar to the 
equilateral model, however, at slightly larger k there is 
considerable support on flattened triangles also. Note 
that for k >, 7.4 fcbump the shape of the bispectrum is 
extremely unusual and is not easily comparable to any 
shape that has been proposed in previous literature. 



V. COSMOLOGICAL PERTURBATION 
THEORY 

In section IIVI we developed an analytical theory of 
particle production and IR cascading during inflation 
which is in very good agreement with nonlinear lattice 
field theory simulations. However, this formalism suf- 
fers from a neglect of metric perturbations and, conse- 
quently, we were unable to rigorously discuss the gauge 







9 As we argued in section III 1 1 the size of the nongaussianity in 
this model is most naturally quantified by evaluating the cum- 
mulants. Here we are interested only in discussing the shape of 
this novel type of nongaussianity. 



FIG. 5: The shape function S(k, kx2, kx-j), defined by (|83[) . 
as a function of the dimensionless quantities X2,X3 which 
parametrize the shape of the triangle. The upper left panel 
corresponds to k = e~ 1 fcb um p, the upper right panel is fc = 
fcbump, the lower left panel is fc = e +1 fcbum P and the lower 
right panel is fc = e +2 fcb um p- In the IR (fc < fcbump) the shape 
of the bispectrum is similar to the equilateral shape, how- 
ever, there is also some support on flattened triangles near 
fc ~ e +1 fcbum P . At larger values of fc the shape is unlike any 
other template proposed in the literature. For illustration we 
have chosen fi 2 — 0. 



invariant curvature perturbation £. Hence, the reader 
may be concerned about gauge ambiguities in our results. 
In this section we address such concerns, showing that 
metric perturbations may be incorporated in a straight- 
forward manner and that their consistent inclusion does 
not change our results in any significant way. We will 
do so by showing explicitly that, with appropriate choice 
of gauge, equations ((55]) and (|3~Tj) for the fluctuations of 
the inflaton and iso-inflaton still hold, to first approxi- 
mation. We will also go beyond our previous analysis by 
explicitly showing that in this same gauge the spectrum 
of the curvature fluctuations, P^, is trivially related to 
the spectrum of inflaton fluctuations, (and similarly 
for the bispectrum). 

To render the analysis tractable we would like to take 
full advantage of the results derived in the last section. 
To do so, we employ the Seery et al. formalism for work- 
ing directly with the field equations 92j and make con- 
siderable use of results derived by Malik in [94|,[95j]. (Note 
that our notations differ somewhat from those employed 
by Malik. The reader is therefore urged to take care in 
comparing our formulae.) 

We expand the inflaton and iso-inflaton fields up to 
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second order in perturbation theory as 



4>{t, x) 



(T,X 
1 

2 



$2X( T , x ) 



(84) 
(85) 



The perturbations are defined to average to zero (5 n (j)) = 
(S nX ) = so that ((f>(t, x)) = <j>{t) and x)) = 0. (The 
condition (x) = is ensured by the fact that m x ^> if 
for nearly the entire duration of inflation.) 

We employ the flat slicing and threading throughout 
this section. With this gauge choice the perturbed metric 
takes the form 



5oo 
9oi 
9ij 



a 2 d, 
a 2 5, 



a 2 (l + 2^i 
Bi 



■to) 



i 

-B 2 
2 



(86) 
(87) 
(88) 



so that spatial hyper-surfaces are flat. Note also that 
in this gauge the field perturbations 6 n <p, S n x coincide 
with the Sasaki- Mukhanov variables [96[ at both first and 
second order. 



rather than by equation (|62|). The power spectrum of 
the gaussian fluctuations is, again, given by (|64]) . The 
correction to the order of the Hankel function v translates 
into a correction to the spectral index: instead of (|6"5|) we 
now have 



1 = 277 — 6e 



(93) 



which is precisely the standard result [1051 ] . 

Thus, as far as the quantum vacuum fluctuations of the 
inflaton are concerned, the only impact of consistently 
including metric perturbations is an 0(e) correction to 
the spectral index n s . 

Let us now turn our attention to the first order fluctu- 
ations of the iso-inflaton. The closed-form Klein-Gordon 
equation for 8\x derived in [94[ can be written as 

S lX " + 2HS 1 x'-V 2 S lX + a 2 [ M 2 + k 4 J 2 (r)] 6 lX = (94) 



This coincides exactly with equation (|31[) , which we have 
already solved. The fact that linear perturbations of x 
do not couple to the metric fluctuations follows from the 
condition (x) = 0. 



Gaussian Perturbations 



B. Nongaussian Perturbations 



In [94[ Malik has derived closed-form evolution equa- 
tions for the field perturbations 5 n (f>, 8 n x at both first 
(n = 1) and second (n — 2) order in perturbation theory. 
Let us first study the gaussian perturbations. The closed- 
form Klein-Gordon equation for 8\<f> derived in can 
be written as 



<5i0" + mSx 







V 2 ^ 



2 2 o 

a m — 3 



Sl<f> 

(89) 



Following our previous analysis we expand the first-order 
perturbation in terms of annihilation/creation operators 
as 



5i(f>(t,x) = 



d 3 k 
(2tt) 3 / 2 



Ok r^— e 

a(r 



h.c. 



(90) 



where h.c. denotes the Hermitian conjugate of the preced- 
ing term and we draw the attention of the reader to the 
the explicit factor of a -1 in our definition of the Fourier 
transform. Working to leading order in slow roll param- 
eters we have 



A; 2 + ^ (-2 + 3?7 - 9e) 



6i<f>k = 



(91) 



This equation coincides exactly with (|61l) and the prop- 
erly normalized solutions again take the form (|63[) . The 
only difference is that the order of the Hankel function, 
u, is now given by 



r? + 3e 



(92) 



Now let us consider now the second order perturbation 
equations. The closed-form Klein-Gordon equation for 
#2</> derived in can be written as 



5 2 <t>" + 2H6 2 cb' - V 2 8 2 < 
= J(r,x) 



2 2 o 

i m — 3 



M„ 



5 2 4> 
(95) 



As usual, the left-hand-side is identical to the first or- 
der equation (|8"9")l while the source term J is constructed 
from a bi-linear combination of the first order quantities 
5\4> and 5%x- In order to solve equation (p?5"f we require 
explicit expressions for the Green function Gk and the 
source term J. The Green function is trivial for the case 
at hand; it is still given by our previous result (I62[) . pro- 
vided one takes into account the fact that the order of 
the Hankel functions v is now given by (|92[) . rather than 
(1621) . In other words, the Green function for the nongaus- 
sian perturbations (|95|) differs from the result obtained 
neglecting metric perturbations only by 0(e) corrections. 

Next, we would like to consider the source term, J, 
appearing in (|95[) . Schematically, we can split the source 
into contributions bi-linear in the gaussian inflaton fluc- 
tuation Si4> and contributions bi-linear in the iso-inflaton 
SiX- 



J — Jrh 



(96) 



The contribution would be present even in the ab- 
sence of the iso-inflaton. These correspond, physically, to 
the usual nongaussian corrections to the inflaton vacuum 
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fluctuations coming from self-interactions. This contri- 
bution to the source is well-studied in the literature and 
is known to contribute negligibly to the bispectrum [9^ |. 
Thus, in what follows, we will ignore J^. 

On the other hand, the contribution J x appearing in 
(|M| depends only on the iso-inflaton fluctuations S\X- 
This contribution can be understood, physically, as gen- 
erating nongaussian inflaton fluctuations Si§ by rescat- 
tering of the produced \ particles off the condensate. 
Hence, the contribution J x may source large nongaus- 
sianity and is most interesting for us. It is straightfor- 
ward to compute JL explicitly for our model using the 
general results of [94|]. We find 

J x = -2 a y(<^0o)(<5ix) 2 

-a 2 (M 2 +.9 2 (0-0o) 2 )(5ix) 2 



+V- 2 (di(5ix)V 2 c>ViX) + V 2 (5ix)V 2 (5!x) 
+S lX 'V 2 6 lX + {Vhx'f ) (97) 



where the upper sign is for <j)' > 0, the lower sign is for 
(f)' < 0. Notice that the contributions to J x on the fourth 
and fifth line of (1^7)) contain the inverse spatial Laplacian 
V~ 2 and are thus nonlocal. These terms all contain at 
least as many gradients as inverse gradients and hence 
the large scale limit is well-defined. In [97[ it was argued 
that these terms nearly always contribute negligibly to 
the curvature perturbation on large scales. 

Let us now examine the structure of the iso-inflaton 
source J x , equation (|97|) . The first line of (|97|) goes 
like a 2 g 2 (<j> — <j>o)(6ix) 2 ■ This coincides exactly with the 
source term in equation (|55[) which was already studied 
in section HVl On the other hand, the terms on the sec- 
ond, third, fourth and fifth lines of (|97| are new. These 
represent corrections to IR cascading which result from 
the consistent inclusion of metric perturbations. We will 
now argue that these "extra" terms are negligible as com- 
pared to the first line. If we denote the energy density in 
gaussian iso-inflaton fluctuations as p x ~ m x (Six) 2 then, 
by inspection, we see that the first line of (p?T|) is para- 
metrically of order p x /\(j} — 4>o\ while the remaining terms 
are or order y // ep x /M p . Hence, we expect the first term 
to dominate for the field values <p = (f>o which are relevant 
for IR cascading. This suggests that the dominant con- 
tribution to J x is the term which we have already taken 
into account in section HVl 

Let us now make this argument more quantitative. 
We assume that p 2 < k 2 , since otherwise particle pro- 
duction effects are exponentially suppressed. Inspec- 
tion reveals that the only "new" contribution to (|9"7} 
which has any chance of competing with the "old" term 
a 2 g 2 {<p— 4>o){5ix) 2 is the one proportional to y/ea 2 g 2 ((j) — 
(f>o) 2 (5ix) 2 /M p (on the second line). This new correc- 
tion has the possibility of becoming significant because 



it grows after particle production, as <p rolls away from 
<f>o- This growth, which reflects the fact that the en- 
ergy density in the x particles increases as they become 
more massive, cannot persist indefinitely. Within a few 
e-foldings of particle production the iso-inflaton source 



term must behave as J v 



, corresponding to the 



volume dilution of non-relativistic particles. Hence, in 
order to justify the analysis of section HVl we must check 
that the term 

Jncw^^aV(0"0o) 2 (<5ix) 2 (98) 

does not dominate over the term which we have already 
considered 



•/old ~«9 (0 - 4>o){8ixY 



(99) 



during the relevant time H At = 0(1) after particle pro- 
duction. It is straightforward to show that 



old 



1 



M p H 1 



1 1 



4>x/e N e N 



(100) 



where N = Ht is the number of e-foldings elapsed from 
particle production to the time when IR cascading has 
completed. Hence, N = 0(1) and we conclude that the 
second, third, fourth and fifth lines of (p?T| are (at least) 
slow roll suppressed as compared to the first line. 

In summary, we have shown that consistent inclusion 
of metric perturbations yields corrections to the inflaton 
fluctuations 8<p which fall into two classes: 

1. Slow- roll suppressed corrections to the inflaton vac- 
uum fluctuations 8i<j) (these amount to changing 
the definition of v in the solution (155)) ). These cor- 
rections have two physical effects. First, they yield 
an 0(e) correction to the spectral index. Second, 
they modify the propagator Gk by an 0(e) correc- 
tion. 

2. Corrections to the source J for the nongaussian in- 
flaton perturbation S2(f>. These corrections are the 
second, third and fourth lines of (|9"7|) ) which, as we 
have seen, are slow roll suppressed. 

It should be clear that neither of these corrections alters 
our previous analysis in any significant way. 



C. Correlators 

So far, we have shown that a consistent inclusion of 
metric perturbations does not significantly alter our pre- 
vious results for the field perturbations. Specifically, 5\X 
is identical to our previous solution of equation (|31[) for 
the iso-inflaton, while 8\4> coincides with the homoge- 
neous solution of equation (|55[) . up to slow- roll correc- 
tions. At second order in perturbation theory, we have 
seen that 

S 2 (j) = J A'G(i-i')J x (i') + 0[(ii(i) 2 ] 
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To leading order in slow roll J x = —2a 2 g 2 ((f> — 4>o)(8ix) 2 
and the first term coincides with our previous result for 
the particular solution of equation ([55]) , The terms of 



order (5i(f) represent nongaussian corrections to the 
vacuum fluctuations from inflation (coming from self- 
interactions of 8<j) and the nonlincarity of gravity). These 
would be present even in the absence of particle produc- 
tion, and are known to have a negligible impact on the 
spectrum and bispectrum [92| . 

We are ultimately interested in the connected n-point 
correlation functions of Scf>. For example, the 2-point 
function ((8(f) 2 ) get a contribution of the form ((81(f) 2 ) 
which gives the usual nearly scale invariant large-scale 
power spectrum from inflation. The cross term (8\4>82(f) 
is of order ((8i(f) A ) and represents a negligible "loop" 
correction to the scale-invariant spectrum from inflation. 
(The cross term does not involve the iso-inflaton since 
8\(j) and 8±x are statistically independent.) Finally, there 
is a contribution ((82(f) 2 ) which involves terms of or- 
der (% 4 ) coming from rescattering and terms of order 
((<5i0) 4 ) which represent (more) loop corrections to the 
scale-invariant spectrum from inflation. Thus, we can 
schematically write 

p (fc) = p; ac (k) [1 + (loops)] + p^(k) 

Here P™ c ~ fc" s_1 is the usual nearly scale invariant 
spectrum from inflation and P™ sc is the bump-like con- 
tribution from rescattering and IR cascading which we 
have studied in the previous section. The "loop" cor- 
rections to P vac (fc) have been stud ied in detail in the 
literature (see, for example, 98l - |i01 |) and are known to 
be negligible in most models. 

We can also make a similar schematic decomposition 
of the bispectrum by considering the structure of the 
3-point correlator ((8(f)) 3 ). Following our previous line 
of reasoning, it is clear that the dominant contribution 
comes from rescattering and is of order (x 6 )- The terms 
involving ((8i(p) 3 ), on the other hand, represent the usual 
nongaussianity generated during single field slow roll in- 
flation and are known to be small 1921. 



D. The Curvature Perturbation 

Ultimately one wishes to compute not the field per- 
turbations S n (f), 8 n x, but rather the gauge invariant cur- 
vature fluctuation, We expand this in perturbation 
theory in the usual manner 



C - Ci + \C2 



(101) 



In [95| Malik has derived expressions for the large scale 
curvature perturbation in terms of the Sasaki-Mukhanov 
variables at both first and second order in perturbation 
theory. We remind the reader that in the flat slicing 
(which we employ) the Sasaki-Mukhanov variable for 
each field simply coincides with the field perturbation 
(i.e. - Q,f, = 8(f) and Q x = 8\)- 



At first order in perturbation theory the iso-inflaton 
does not contribute to the curvature perturbation (since 
(x) = 0) and we have 



(102) 



At second order in perturbation theory the expression 
for the curvature perturbation is more involved. Using 
the results of (95j | and working to leading order in slow 
roll parameters we find 10 



u 



8 2 (j> ■ 



8j£ 

m 



(103) 
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+ ^2 [(Six') 2 + a 2 (^+9 2 v 2 t 2 (r))(8 1 xf 



+ ^[^r+a 2 m 2 (8^) 2 ] 

Let us discuss the various contributions to this equation. 
The third line contributes to the nongaussianity of the 
vacuum fluctuations durin g in flation. These terms are 
known to be negligible [5m58l l92j | and, indeed, one may 
explicitly verify that (|103p would predict Jnl ~ 0(e,7]) 
in the absence of particle production. 

Next, we consider the second line of (|103[) . This repre- 
sents the direct contribution of the gaussian fluctuations 
81X to the curvature perturbation. This contribution is 
tiny since the x particles are extremely massive for nearly 
the entire duration of inflation and hence Six ~ a -3 / 2 
(see also Q for a related discussion). The smallness of 
this contribution to £ can be understood physically by 
noting that the super-horizon iso-curvature fluctuations 
in our model are negligible. 

Finally, let us consider the contribution on the first 
line of (I103p . This contribution is the most interesting. 
To make contact with observations we must compute the 
curvature perturbation at late times and on large scales. 
In section HVl we have already shown that S n (j) is constant 
on large scales and at late times for both n = 1 and n — 2. 
This is the expected result: the curvature fluctuations 
are frozen far outside the h orizon and in the absence of 
entropy perturbations [lOJ]. 11 Hence 82$ is completely 
negligible and the first term on the first line of (|103l) must 
dominate over the second term. We conclude that ,at 



We have dropped a sp uriou s additive 2C^ which stems from using 
the Malik and Wands 102] definition of the curvature perturba- 
tion, rather than the definition employed by Lyth and Rodriguez 
|103(1 and also by Maldacena [57I . (See also fiojl . *) 
Note that, in some cases, the curvature fluctuatio ns m ay evolve 
significantly after horizon exit [106| . I107H . (See also |10S( . ) This is 
a concern in models where there are significant violations of slow- 
roll. In [l[ we have already shown that the transient violation of 
slow roll has a negligible effect on the curvature fluctuations in 
our model; see also [72j . Hence, the result f n ~ S n <j> ~ const far 
outside the horizon is consistent with previous studies. 
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late times and on large scales, the second order curvature 
perturbation is very well approximated by 



<2 



u 



-5 2 4> ■ 



(104) 



In summary, we have shown that the power spectrum 
of curvature fluctuations from inflation in the model @ 
is trivially related to the power spectrum of inflaton fluc- 
tuations 



p c (k) - 



i 



2eM 2 



(105) 



at both first and second order in cosmological perturba- 
tion theory. This relation is valid at late times and for 
scales far outside the horizon. The curvature spectrum 
(|105[) may be written as 



P ( (k) = P^ ac (k) [1 + (loops)] + P c resc (£;) 



(106) 



The power spectrum of the inflaton vacuum fluctuations 
agre es with the usual result obtained in linear theory 



r>vac 



(k) 



H 2 



8n 2 eM 2 



aH J 



(107) 



In (|106[) we have schematically labeled the corrections 
arising from the third line of (|103|) and the source as 
"loop" . These are nongaussian corrections to the inflaton 
vacuum fluctuations arising from self-interactions of the 
inflaton and also the nonlinearity of gravity. Such correc- 
tions are negligible. The most interesting contribution to 
the power spectrum (|106|) is due to rescattering, P™ sc (k). 
This quantity is proportional to our previous result (|77[) . 

In passing, notice that the bispectrum (defined by 
(|80p) of inflaton fluctuations will differ from the bispec- 
trum B of the curvature fluctuations (defined by (jYTJ) ) 
only by a simple re-scaling: 



B(h 



H 



1 



(2e) 3 / 2 M p 3 



B^{h) (108) 



The dominant contribution to B^ comes from rescatter- 
ing effects and scales as {82^) ~ (<5iX 6 )- 

The analysis of this section justifies our neglect of met- 
ric fluctuations in section HVl 



VI. CONCLUSIONS 

In the context of a realistic microscopic frame-work, 
we might generically expect the inflaton to couple to a 
large number of fields whose energy density does not play 
any important role in driving inflation. Such couplings 
can lead to isolated bursts of particle production during 
inflation. The associated observational signatures pro- 
vide a rare opportunity to learn about how <f> couples to 



other species, as opposed to the self-coupling information 
which is encoded in V(<p). In this paper we have consid- 
ered a simple example of this effect which is dynamically 
rich and derivable from realistic particle physics models, 
such as string theory. 

Inflationary particle production leads to features in the 
primordial curvature fluctuations via the mechanism of 
IR cascading. This process is interesting in its own right: 
it is qualitatively different from other mechanisms in the 
literature (in that we do not rely on the quantum vacuum 
fluctuations of some light iso-curvaturc fields) and the 
underlying dynamics are relevant for preheating, mod- 
uli trapping and non-equilibrium QFT more generally. 
Moreover, particle production and IR cascading lead to 
a variety of novel observable signatures, including local- 
ized features in both the spectrum and bispectrum of the 
cosmological fluctuations. 

In this paper we have extended previous work [3, Hj] on 
inflationary particle production in two directions. Firstly, 
we have developed an analytical theory of particle pro- 
duction and IR cascading during inflation, which is in 
excellent agreement with lattice field theory simulations. 
This formalism helps to clarify the underlying physics of 
the mechanism, and provides a crucial cross-check on our 
numerical methods. 

Our second main result has been a more detailed in- 
vestigation of the nongaussian signature associated with 
particle production and IR cascading. The bispectrum 
in this model is rather unusual: it peaks only for trian- 
gles with a size comparable to some characteristic scale. 
We have argued that the magnitude of this type of non- 
gaussianity is best characterized by studying the mo- 
ments of the PDF. For realistic values of the coupling, 
the skewness of the PDF is quite large. For example, 
with g 2 ~ 0.01 the power spectrum for our model is com- 
patible with all observational data [2| while the skewness 
of the PDF is equivalent to what would be produced in a 
local model with /Jv£ lv ~ — 53. This value is somewhat 
larger than current observational bounds, suggesting that 
nongaussianity from inflationary particle production may 
be observable in future missions. However, we stress that 
the nongaussian signature in our model is quite differ- 
ent from what would be expected for a local model with 
C - C 5 + lf N t V [Q - (Cg)}- ^ particular, the higher 
order cummulants (such as the kurtosis) are different, as 
are the shape and running of the bispectrum. 

Note that, if it were to be detected, the nongaussian 
signature from IR cascading must be correlated with an 
observable feature in the power spectrum and also with 
signatures in polarization. Hence, it should be possi- 
ble to robustly rule out the possibility that massive iso- 
curvature particles were produced at some point during 
the observable range of e-foldings of inflation. 

The nongaussian signature predicted by inflationary 
particle production is rather complicated as compared 
to the local or equilateral models. However, the under- 
lying field theory description of our model is extremely 
simple and rather generic from the low-energy perspec- 
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tive. In order to obtain large nongaussianity it was not 
necessary to fine-tune the inflaton trajectory or appeal 
to re-summation of an infinite series of high dimension 
operators. Indeed, the only "tuning" which is required 
for our signal to be observable is the requirement that 
<fi = <f>0 during the observable range of e-foldings. We 
believe that this type of nongaussianity is very natural 
and merits further investigation from the observational 
perspective. 

There are a variety of directions for future studies. 
From the theoretical perspective, it would be interesting 
to explicitly generalize our results to more complicated 
models with particle production during inflation (such as 
SUSY models, higher spin iso-inflatons and phase tran- 
sitions). There are also a wide range of interesting phc- 
nomenological possibilities. Varying the location of the 
feature we can have a variety of possible signatures for 
the CMB and LSS. We expect that IR cascading will 
also have implications for the spectrum of gravity waves 
from inflation and also primordial black holes. We could 
imagine superposing multiple bursts of particle produc- 
tion to obtain an even richer variety of signatures. It 
would be interesting to construct a simple, separable es- 
timator for the bispectrum from IR cascading which can 
be confronted with observational data in order to obtain 
explicit constraints on the underlying model parameters. 
We leave these possibilities for future investigation. 
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APPENDIX A: Detailed Computation of P(k) 



In this appendix we discuss some of the technical de- 
tails associated with the computation of the renormalized 
power spectrum (|77|) . First, notice that using ([50| and 
(|51l) we can write the quantity appearing in each renor- 
malized Wick contraction as 



X k (r)xUr r )-Mr)r k (r') = 
1 I 1 



kl y/a(r)a(r') 



n k cos 



klt^r) klt\r') 



n k sin 



k^(r) kit\r>) 



(A-1) 



This work is dedicated to the memory of L. Kofman, 
who played a key role in initiating this line of research. 



where the occupation number n k is defined by (149|) . Plug- 
ging 

CSH into dTTJ) we find 



g 2 k 3 
8tt 5 



d 3 k'n k _ k in k i x / dr' dr 
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kie{r<) , klt\r")- 
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x / dr'dr' 
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a(r) a(r) 



rfc^ 2 (r') klt\r") 



k^(r') , klt\r") 



(A-2) 



Notice that the time and phase space integrations A.l. Time Integrals 

in (|A-2[) decouple. This is the key simplification 
which makes an analytical evaluation of this expression 
tractable. Let us consider these integrations separately. 



All of the integrals over conformal time that appear 
in (|A-2I) can be expressed in terms of two characteristic 
integrals which we call I\ and li- Explicitly, these are 
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defined as 

h(k,T) 



a(r) 



dr'G k {T-T')e lk " t2 ^ (A-3) 



h{k,T) = t~t" ( dr'G k {T-T') 
a{T) J 



(A-4) 



The second characteristic integral, I2, can be evaluated 
analytically. However, the resulting expression is not par- 
ticularly enlightening. Evaluation of I\, on the other 
hand, requires numerical methods. 

Let us now show how the various integrals appearing 
in (|A-2j) may be re- written in terms of ii, 12- First, 
consider the first line of (|A-2I) where the following integral 
appears: 

, , , „ G k {T-T')G k (T-T") 

drdr — 

a(r) a(r) 

k 2 J 2 {r>) klt*{r") 



cos 



\h(k, T )\ 2 , HKrf 



Next, consider the second line of (|A-2[) where the follow- 
ing integral appears: 

, , . „ G k (T - T 1 ) G k {T - T") 

drdr — -— 

a(r) o(r) 

k 2 t 2 (r') , k 2 t 2 (r") 



2 2 
Re[h(k, T ) 2 ) I 2 (k, T ) 2 



(A- 



2 2 

Finally, consider the fourth line of (|A-2[) where the fol- 
lowing integral appears: 

dr'dr^-^^r:' 1 

o(r) a(r) 
fc 2 i 2 (r') fc 2 t 2 (r")l „ ; „ rfc 2 i 2 (r') , fc 2 t 2 (r") 



2 2 

Im [/!(*, T)J 3 (fc,T)] 



This formula is valid when a, 6 are positive real num- 
bers. Notice that this expression is symmetric under in- 
terchange of a and b. 

The phase space integral in the first line of (| A-2|) is 
computed by a trivial application of the identity (|A-8|) : 



d 3 k'n k - k >n k > = 



k 3 



2y/2 



(A-9) 



However, the remaining phase space integrals appearing 
in (| A-2|) cannot be obtained exactly in closed form be- 
cause they contain terms like yl + n k > where the gaus- 
sian factors appear under the square root. In order to 
deal with such expressions, we note because n k <C 1 over 
most of the domain of integration, it is reasonable to re- 
place yl + n k i = 1 + n k > /2. Let us now proceed in this 
manner. The phase space integral on the second line of 
(|A^2|) is: 



d k'^n k - k :n k > y/1 + nk-vv'l + n k , 
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Finally, consider the phase space integral on the third 
line of {OJ: 



d 3 fc' 
d 3 /fc' 
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In the expressions (|A-6|) and (|A-7[) the notations Re and 
Im denote the real and imaginary parts, respectively. 



A. 2. Phase Space Integrals 

As a warm-up to the subsequent calculation consider 
the following integral: 



d k TLfe fe/7lfc/ 

J d 3 k' exp [-a7r|k-k'| 2 /fc 2 ] exp [-6vr|k'| 2 /A: 2 ] 
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2V2 



e 3^ 2 /(2fe^) exp 



e ^/(2kl) exp 



irk 2 
'3fc 2 

37Tfc S 
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We have verified the formulae (|A-10IA-11[) numerically. 
In both cases that the numerical results agree with these 
semi-analytical expressions up to the percent level. 

We can now, finally, insert the results (|A-5IA-6IA-7[) 
and (|A-9IA-10IA-11|) into the expression (|A-2|) . Doing so, 
we arrive at our main analytical result, which is equation 

ca. 
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APPENDIX B: Evaluation of the Bispectrum 

In this appendix we discuss some of the technical 
details associated with the explicit evaluation of the 
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renormalized bispectrum (|77[) . To render the analysis 
tractable we will work in the flat-space limit H — > 0, 
which is sensible since the process of IR cascading takes 
only a single e-folding and, moreover, this approximation 
was shown to yield sensible results for the power spec- 
trum in [l[ . Neglecting the expansion of the universe we 
have t = t and the retarded Green function (|66[) becomes 



Gk{t - if) = Q{t n sin [n k (t - f )] 

ilk 



(B-l) 



with Q(x) the Heaviside function and ilk = Vk 2 + m 2 . 
In this limit the characteristic integrals I±(k,t), l2(k,t) 
defined by ([TSijl and (f79j) can be computed analytically. 
For Ii(k,t) we find 
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while, for ^(fe,/;), we have 



/ 2 (M) 
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(Note that our definition of Ji , J2 differs from [H by a fac- 
tor of ^fc 1 -) Finally, the renormalized Wick contraction 
(|A-1[) also simplifies in the limit H — > 0: 
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where the occupation number n k is defined by (fTU 



Inserting (|B-1|1 and ()B-5j) into (81]) we find the follow- 
ing expression for the renormalized 3-point correlation 
function: 



4 5 3 



(2 , )9/ ^ (3) (ki + k 2 + k 3 ) n 7^ s - ^ c* - h 



x / d 3 p 



n kl ~ p cos 



(Mi) 2 (M 2 ) ; 



-n Sill , 

p v 2 



n fc3+p cos 



(M2) 2 (M 3 ) : 



n n cos 



") +"fe(-pV 1 + n fci- 



(Mi) 2 (Ma) s 



(M 2 ) 2 (Ms) 1 



2 2 

2 (h + \2 



(Mi) 2 (Kt 3 ) 2 \ 1/2 n ——. ( (Kh) 2 (M 3 ) 



(B-6) 



It only remains to expand out the expression (|B-6|) and 
evaluate the various integral which arise. As in the case of 
the 2-point function, the phase space and time integrals 
decouple, making an analytical evaluation tractable. Let 
us consider the various integrals which arise separately. 



B.l. Phase Space Integrals 

First, let us introduce a notation for the fundamental 
phase space integral which arises 



K a ,b,c = J d 3 p n a kl _ p n\ 3+p n c p 
kl 



■ exp 
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x exp 
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where we have used the fact that k x + k 2 + k 3 = 0. To 
evaluate integrals containing radicals suc h as yj \ + n p we 
use the same trick as was employed for (|A-10I) . That is, 
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we approximate: 

J d 3 p n kl ^ p n k3+p np /2 y/l + n p 
= J d 3 p n kl ^ p n k3+p n p /2 + i J d 3 p n kl 
= #1,1,1/2 + 9 #1,1,3/2 



that this gives a good approximation to the exact result. 



B.2. Time Integrals 



The evaluation of the time integrals appearing in (|B-6[) 
(B-8) is a straightforward generalization of the results pre- 
sented in appendix A. Let us introduce some notations 
and similarly for the other combinations which arise in for the various combinations of the fundamental integrals 
the expansion of (|B-6[) . We have checked numerically I\ and I 2 that will appear in the final result: 
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B.3. The Full Bispectrum 

We are now finally in a position to write out an ex- 
plicit expression for the renormalized 3-point function 
of the inflaton fluctuations generated by IR cascading. 



That expression is given below, in terms of the various 
that were defined explicitly in equations (|B-7|) and (IB-9p ~ 
(|B-16j) . As promised, the explicit result for the 3-point 
correlation function is cumbersome and not entirely en- 
lightening. 
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